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Abstract 

The motion of a continuum of matter subject to gravitational interaction is classically de- 
t^- ■ scribed by the Euler-Poisson system. Prescribing the density of matter at initial and final times, 

we are able to obtain weak solutions for this equation by minimizing the action of the Lagrangian 
Qh I which is a convex functional. Through this variational formulation, the reconstruction problem 

■ becomes very similar to an optimal transportation problem. Then we see that such minimiz- 

ing solutions are consistent with smooth solutions of the Euler-Poisson system and enjoy some 
' special regularity properties. 

_ 1 Introduction 

> 

O . The Euler-Poisson system describes the motion of a self-gravitating fluid. It is used in cosmology, 
^ | to model the evolution of the primitive universe. In the classical (non-relativistic) description, the 
^0 gravitational field generated by a continuum of matter with density p, is the gradient of a potential 
^ ■ p satisfying the Poisson equation. The system is thus the following 

§: d t p + v-(pv) = o, 

d t (pv) + V • (pv <g> v) = -pVp, 
Ap = p, 



X 

5-H 



where v is the velocity field. The model is expected to be valid at scales where no collisional effects 
enter into account, and no multi-streaming occurs, otherwise the Vlasov- Poisson system would offer 
?3 \ a preferable description. 

Assuming that the Universe is described by this set of equations, we are now interested in solving 
the so-called reconstruction problem, that we introduce here. We can infer, from red-shift catalogs, 
the repartition of matter in the present Universe. We also know that at times very close to the Big 
Bang, the matter was highly concentrated, but with very small relative fluctuations of density. From 
this knowledge, can we reconstruct the intermediate states of the Universe, as well as the initial and 
present velocities ? The problem amounts to reconstruct a solution of the above Euler-Poisson system 
based on the knowledge of the initial and final density fields, and seems ill-posed. Another condition, 
known as slaving reduces the number of unknowns, and renders the problem well-posed: we impose 
that the reconstructed solution has a potential velocity field, i.e. v = V0 for some potential <p. This 
paper proposes a way to solve the reconstruction problem, by transforming it into a minimization 



1 Laboratoire J.A.Dieudonne, Universite de Nice-Sophia- Antipolis & Ecole polytechnique federale de Lausanne 



1 



problem, in many ways similar to an optimal transportation problem. Most of the results will be 
obtained by use of the Monge-Kantorovitch duality, a tool widely used in optimal transportation 
problems. 



1.1 Equations of motion in co- moving coordinates 

The system of equations we will look at, is not exactly the one displayed above. We will first 
express the equations of motion taking into account the global expanding motion of the Universe. 
We consider the expansion factor a(t), a scalar function of time. Following the global expansion, 
particles positions x(t) are given by x(t) = a(i)xo and xo is called the co-moving coordinate. The 
local velocity corresponding to this uniform expansion is then v u (t,x) = (a/a)x. The mean value 
of the density p m is supposed to be close to the critical density, defined as the highest value of p m 
allowing an infinite expansion; this condition is known to imply that a(t) = (t/t ) 2 ^ 3 . Performing 
then the change of variables 

P — PPm, v — (d/a)x + a(t)v, 
P = p m +P, r = (t/t ) 2/3 , 

we obtain the new Euler-Poisson system 

d T p + V • (pv) = 0, 

(1) d T v + v- Vu = -(3/2r)0 + Vp), 

Ap = (p - l)/r. 

Note that the homogeneity condition at time reads p(0) = 1. The slaving constraint appears as 
a necessary condition for the right hand side of the momentum equation not to be singular as time 
goes to 0. In this work, we will consider a simplified version of this system where we do not include 
the time dependence in the coupling between p and p, as well as the drag term — j-v in the second 
equation. We will also restrict ourselves to the flat torus T d = IR d /Z d with p m = 1. In this framework 
the Euler-Poisson system hereafter referred to as (E — P) takes the following form: 

(2) d t p + V-(pv) = 0, 

(3) d t (pv) + V • (pv ® v) — -pVp, 

(4) Ap = p - 1, 

with the additional constraint 

p(.,x)dx = 1. 

■d 

The only modification from the first set of equations is the neutralizing background effect, that 
transforms the Poisson equation Ap = p into Ap = p — 1. One of the main consequences is that in 
this model, p = l,t> = 0isa solution. This corresponds in the physical coordinates to a uniform 
expansion. As explained in paragraph 11.61 below, the results and techniques used for this simplified 
system extend naturally to the full system (|TJ) . Note that this form is the cosmological one, i.e. that 
the potential is attractive. In the case of a repulsive potential (used for the description of a plasma) 
the associated Poisson equation would be Ap = — [p — p m }. 
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1.2 Definition of the reconstruction problem 



Given the (E — P) system, one can try to solve the Cauchy problem, i. e. given p and v at time t = 
find a solution to (0EHH) oia a. time interval [0, T[. Another approach is to look for a solution over 
the time interval [0, T] satisfying the two conditions: 

(5) p\ t =o = Po, 

(6) p\ t =T = Pt- 

This approach has been used by Brenier in j^J, [7] for the incompressible Euler equation and allows 
to introduce variational techniques. Indeed the system (0 El EJ) is hamiltonian, with Hamiltonian (or 
energy) given by: 

H{p, v) = I I \v(t, x)\ 2 dp(t, x) - \Vp{t, x)\ 2 dx. 
i jjd 

Solutions of hamiltonian systems are critical points for the action of the Lagrangian, here defined by 

(7) I(p,v,p) = ^ [ I \v(t,x)\ 2 dp(t,x) + \S7p(t,x)\ 2 dxdt, 

* JO JT d 

under the constraints (J21 EJ EJ El) • Expressed in terms of (p, J = pv) the action becomes 



1 - ' ' - 

F(p,J,p) = - 



J7 d 



dp(t,x) + \ Vp(t,x)\ dxdt. 



We will see that this function is convex in (p, J, p), therefore the critical point will necessarily be a 
minimizer. The goal of this paper will thus be to solve the following problem: 



Problem 1.1 Find p,v,p such that 

I(P'V,p) = MI(p,v,p) 

over all the triple (p, v,p) satisfying 

d t p + V- (pv) = 0, 
Ap = p - 1, 
p\t=o = Po, 
p\t=T = Pt- 

The problem is here formulated in a very vague way: we do not mention in what space lie p and v 
when we perform the minimization. This will be made precise in the next subsection. 



1.3 Motivations 

The interest of studying this boundary problem is twofold. The natural motivation is its direct 
application in cosmology, for the reconstruction of the early Universe. A first approach had been made 
in jH], where the authors solved a variant of the problem, assuming the Zel'dovich approximation. 
This approximation turned the reconstruction problem into an optimal transportation problem with 
quadratic cost. The dual Monge-Kantorovitch problem was solved numerically using an algorithm 
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due to Bertseakas. A more detailed discussion about the physical aspects of the problem was also 
presented in that paper (see also |16j). The reader can also refer to the PhD thesis of J. Bee ([2], 



On the other hand, from a mathematical point of view, it has been observed first by Brenier 
in [Zj in the case of the Euler incompressible equation and by Evans and Gomes in ^3] in the 
smooth finite dimensional case that solutions of Hamiltonian flows that minimize the action of the 
Lagrangian present additional interesting regularity features. Some aspects of the present work can 
be seen as a continuation of their contribution in one special case of infinite dimensional Hamiltonian 
system. In a very vague way, one of the common features of between those three approaches is the 
following: the dual formulation of the problem lets appear a Hamilton- Jacobi equation, coupled with 
a transport equation. Solutions of Hamilton- Jacobi equations are not expected generally to present 
much regularity (usually they are not differentiable), due to the formation in finite time of caustics. 
However, thanks to our variational approach, the solutions found here will behave in a nicer way 
than one can expect. On a sets that contains the dynamics (i.e. that contains the support of p) 
the solution of the Hamilton- Jacobi equation will be differentiable, and even C l ' LXogL (see point 7 
of Theorem 12.3)1 . Moreover, the dynamics will be reversible on this set. Therefore, this study shows 
some links between the weak KAM theory developed by Fathi ([15]), Evans (previous reference and 
[T3*j). optimal transportation and Monge-Kantorovitch duality. 

Another interest of this study is to generalize the "Lagrangian minimizing" approach developed by 
Benamou and Brenier in [H], in which the authors gave a continuum mechanics interpretation of the 
Monge-Kantorovitch problem involving the concept of interpolation between two measures, induced 
by the Wasserstein distance. This interpolation was also introduced earlier by McCann in |2Tj to 
develop the useful concept of displacement convexity, that will also appear in this case. Otto in [22] 
also used it to endow the set of probability measures with a formal Riemannian metric, in which 
the interpolation plays the role of geodesies, then allowing rich interpretations of some dissipative 
equations in terms of gradient flows. Here our variational problem induces an interpolation that has 
somehow more regularity than the one of [2] where the Lagrangian is only 



Indeed we will see that the additional Dirichlet term forces the intermediate densities to be in L°°{T d ) 
independently of the initial and final densities. One can give the following heuristic interpretation of 
this result: the gravitational force being attractive, the system has to go first through an expansion 
so that a concentration does not appear before the final time T; this fact is expressed through the 
differential inequality ()4*3|) . Meanwhile, some interesting displacement convexity properties, similar 
to the one found in [21], will also appear. 

The techniques we will present here can be naturally adapted to minimize the functionals 



provided T is a convex functional in p. Examples are T(p) = p 7 , yielding solutions of the gas 
dynamics system, with attractive pressure term. This system has an application in the theory of 
large deviations for random matrices, see [T8j . 

Another example is the simple case ^(p) = pV, for some potential V(t,x). This yields solutions 
of d t (pv) + V ■ (pv ® v) = -pW. 

We finally mention that the first steps of this study had been already done in E. Camalet's PhD 
lOj under supervision of Y.Brenier. 



p.11,12). 
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1.4 Organization of the paper 



The paper is organized as follows: we first give a slightly different form to the minimization problem 
( Problem I1.3J1 . which we show is equivalent to the formulation of Problem 11.11 We then state our 
results: existence, uniqueness of the minimizer, optimality equations (Theorem 12 .2|) . and regularity 
properties of the optimal path (Theorem 12. 3j) . 

The rest of the paper is devoted to the proof of those results: in section 3 we show existence 
and uniqueness, in section 4 we derive the important formula (|20p. this will yield that the optimal 
path solves the Euler-Poisson system (section 5), and some partial regularity results (section 6). In 
section 7, we then investigate the consistency of our minimizing solution with other solutions of the 
Euler-Poisson system. This ends the proof of Theorem 12.21 The section 8 is devoted to the proof 
of Theorem 12.31 for this we introduce a time discretization of the problem (Problem 18. lj) . For this 
approximate problem we are able to obtain rigorous estimates, that yield the result when the time 
step goes to 0. 

Formal proofs are given that present all the arguments necessary for the main results: in paragraph 
4.1 for Theorem 12 .2\ and in paragraph 8.1 for Theorem 12.31 



1.5 Precise definition of Problem ll.il 

We introduce the domain D = [0, T] x T d . We also define the flux of matter J by J = pv. Given 
J G M. d , p G M + we use the fact that 

( +oo if p = 0, J^O, 

(8) sup {pc + J ■ m} = < if J = 0, 

cm,m£R d ,c+\m\ 2 /2<0 ] j2 j^p if p > 0. 

Notice that as a supremum of affine functions, this is a (possibly infinite) convex functional in (p, J). 
Given po,px as in Theorem 12.21 the functional I can thus be formulated as 

(9) I(p,J,p) = sup < / c(t, x)dp(t, x) + m(t, x) ■ dJ(t, x) > + - / |Vp(t, x)\ 2 dtdx, 

c+|m| 2 /2<0 I Jd ) 2 J D 

where the supremum is taken over all (c, m) G C{D) x (C(D)) d . As we shall see in the next 
proposition, this formulation is consistent with the formulation ((7|) in the case where v G L 2 (D,dp) 
and well defined (although leading to possibly infinite value) for p G C([0, T]; V(T d ) — w*), J G 
(Ai(D)) d , Vp G L 2 (D) where Ai(D) denotes the set of bounded measures on D and V(T d ) the set 
of probability measures on T d . 

Proposition 1.2 Let p be a measure [0, T] x T d , v be a dp measurable vector field, and p be a 
measurable function on [0, T] x T d . Let the junctionals I and I be defined by (0) and (0) re- 
spectively, with the convention that I = +oo whenever v ^ L 2 (D,dp) or Vp ^ L 2 (D,dx). Then 
I(p,v,p) = I(p,v,p) as junctionals valued in } — oo, +oo]. 

Proof. First suppose that J is not absolutely continuous with respect to p. Then there exists a 
set S such that | <7| (aS) > and \p\(S) = 0. Then one constructs a sequence c e ,m t such that 

J c e dp + m e • d J — > +oo. 



5 



Indeed there exists a sequence O e of open sets containing S such that |p|(O e ) < e and |J|(O e ) > 5 
with 5 > fixed. Then for each e, there exists f e E {C™{O e )) d , ||/ e || L °° < 1 such that J f e -dJ> 5/2. 
We then have, taking m e = A/ £ , c e = — \m e \ 2 /2, 



J -^\\f t \ 2 dp + A/ e • d.J > X5/2 - \\/2 



which is greater than 5e 1 ^ 2 /2 — 1/2 for A = e 1 ^ 2 . Letting e go to 0, we conclude that definitions 
(0) and © yield +oo. 

v 

If J << p, J = pv for some p measurable function v; we consider vp = — max{\v\,P}. By 
standard smoothing arguments, there exists a smooth sequence (c n = —\m n \ 2 /2,m n ) such that 

J c n dp + m n ■ vpdp = J dp\vp\ 2 /2 — J dp\m n —vp\ 2 /2 

dp\v P \ 2 /2, 

/i ^ 1 2 
dp > 0, and this implies 

J cpdp + ra P ■ vdp = J dp\vp\ 2 /2 — / dp\v p — m P \ 2 /2 + J dp(v — vp) ■ v p 

> J dpv-v p - dp\v P \ 2 /2 - e(P) 

> [dp\v P \ 2 /2-e(P), 



with e(P) -> as P -> +oo. Hence J > L 

Then, if J p|f | 2 is bounded, for any c < — \m\ 2 /2, 



112 I 12 112 

\v\ , v — m\ / , f 



cdp + m ■ vdp < I dp— dp < / dp—^-. 



This shows I < I whenever / < +oo, and ends the proof of Proposition II. 21 

The new formulation of the Problem 11.11 is then: 
Problem 1.3 Minimize 

I(p,J,p) = sup < / c(t, x)dp(t, x) + m(t, x) ■ dJ(t, x) 

c+|m| 2 /2<0 yj D 



□ 



+- / \Vp(t,x)\ 2 dtdx 
2 Jd 



among all (p,J,p) that satisfy p G M.(D) H C([0, T]; V(T d ) — w*), J E (M(D)) d , VpEL 2 (D), and 
satisfy in the distribution sense 

(10) d t p + V ■ J = 0, 

(11) Ap = p-1, 

(12) pit = 0) = po, 

(13) pit = T) = p T . 
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We denote 



K=MI(p,J,p) 



p,J,p 



among all such (p, J,p). 

1.6 Variational problem for the full equations in co-moving coordinates 

As we have seen above, the full system for which we want to solve the reconstruction problem is the 
following 



The same techniques as the one that will be exposed in this paper adapt to this minimization problem. 

2 Results 

2.1 Notation 

The space dimension will be denoted d. 

We shall hereafter use the notation A: B = trace (A t B) for A, B two d x d matrices. 
For a function (p : R d — > R m we denote by Dip its first derivative matrix equal to ( 1^ ) 1 < i < . 



( d T p + V-( P v) = 0, 

< d T v + v- Vv = -3/(2r)(u + Vp), 

[ Ap = (p - l)/r. 



Solutions for this system can be sought as minimizers of the action 




under the constraints 



d T p + V • (pv) = 0, 

p(r = 0) = 1, p(r = T) = p T . 




2.2 Definition of weak solutions for (E — P) 

Definition 2.1 A triple (p,v,p) is said to be a weak solution of (E — P) if: 

1. P E L 2 ([0,T];//- 1 (T <i ))nC([0,T];7 ? (T <i ) - w*), v E L 2 (D,dp), 

2. for any y? = (p j ) m ..d] G (C c °°(]0, T[xT d )) d one has 



(14) 




d t ip ■ v dp + Dip : v <g> v dp — p ■ Vp 
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3. for any <p G C°°([0,T] x T a 



d t (fdp + V(f ■ vdp = / p T (f\t=T- / Po<p|t=o, 

[0,T]xT d JT d JT d 



I {dp - l)tp + Vp • V<p = 0. 

</fO,TlxT d 



'[0,T]xT d 

Equation (jTj| is equivalent to equation Q for smooth p using the identity 

(1 + Ap)Vp = V ■ (Vp <8> Vp) - ^V|Vp| 2 + Vp. 

The right hand side of this identity is well defined in the sense of distribution if we only know that 
Vp G L 2 (D). 

2.3 Statement of the Theorems 

We first have an existence/uniqueness result for the minimizer of the action I defined in (|9j). Note 
that from Proposition II .21 this minimizer yields also the minimizer of the action / defined by (J7J). 

2d , 

Theorem 2.2 Let po,Pr be two probability measures in L d + 2 (T ), then there exists a unique 
(p, J,p) G (M(D) x (M(D)) d x L 2 ([0, T\; H x {T d ))) with Ap = p - 1 in V minimizer of the Problem 
\1.3\ The flux J has a density v with respect to p, (p,v,p) is a weak solution of the Euler Poisson 
system (E — P) in the sense of Definition \2.1\ and coincides with any smooth solution of {E — P) 
satisfying M 6 A and having a potential velocity; such solution is therefore unique. Moreover 

1. there exists <f> G Lf oc (]0,T[; H l (T d ))nL™ c {]0,T[xT d ) such that v = V<f) dp a.e. and we can thus 
extend the definition of v to all ofT d as a function belonging to Lf oc (]0,T[xT d )) , 

2. any such extension satisfies 

(15) / / \v(t,x + y) - v(t,x)\ 2 dp(t,x) < C T \y\ 2 



!jd 

for all t in ]0,T/2], y in R d , 
3. the density p belongs to Lf oc {]0, T[xT d ) n C(]0, T[; LP) for any p G [1, 3/2[. 
Then we have the regularity result: 

2d . _ 

Theorem 2.3 If po and px are in L d + 2 the unique solution (p, J, p) of Problem M .31 has the following 
regularity properties: 

1. The density p is in L% c (}0, T[xT d ) n CQO, T[; L k {T d )) for every 1 < k < oo: for every 
t g]0, T/2[ there exists C T such that for every t in [r, T — t], 

IIp(^ -)IU°°(T d ) < C T , 

and there exists C such that 

-C(l + \)<J t (\\p(t, -)\\ Lkm ) < C(l + ^), 
moreover the constants C T , C are independent of the choice of po and pt- 
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2. The velocity v = V0 can be chosen in Lf£ c (]0, T[xT d ), this bound is also independent of the 
choice of p and px- 

S. nejunction, k>l, [ [p^Mt^dx are convex vith respect to time. 

JT d Jfd 

4- The velocity potential <p can be chosen in W lt £°(\Q, T[xT d ) and to be viscosity solution of 
d t (p+ ||V0| 2 +p = on every [s,t] c]0,T[. 

5. If pt is in L p (T d ) with p > d then point 4 extends up to t = T . 

6. One can also choose (ft such that {jp,q){t,x) = (— 4>,p)(T — t,x) is a viscosity solution of 
dtip + ^\Vip\ 2 + 1 = on every [s,t] c]0, T[ and point 5 applies. Moreover both choices will 
coincide dp almost everywhere. 

7. For each t e]0, T[, there exists a closed set S t of full measure for p(t) such that 4>{t) is differ- 
entiate with respect to the space variable at every point of St- Moreover, for all t G [r, T — r], 
r > 0, for all (x,y) in S t , with \x — y\ < 1/2, 

|V0(t, x) - V0(t, V)\ < C(r)\x - y\ l g(— !— A 

\x - y\ 

The reader can refer to the books of Evans ^2] and Barles pQ for the definition of viscosity solution 
of d t <p+\\V<P\ 2 +p = $. 

Remark 1. The assumption that the final and initial densities are in L~ d ~+5{Y d ) is technical: it 
allows us to show that there exists at least one admissible path with finite action transporting p on 
pr, see section 13.0.11 Actually all the results are true assuming that there exists a (p, J,p) satisfying 
all the constraints fp~0| . . . fT3|) such that I(p, J,p) is finite. 

Remark 2. The bound (|T5|) is a finite difference version of the formal (but non rigorous since p has 
no regularity ) assertion J[ rT _ T ] xX d dp\ Vf | 2 < +oo. See jl] where the authors look at an appropriate 
definition of the tangent space related to a measure. 

Remark 3. One may observe that the regularity obtained for p is stronger in the second theorem 
than in the first. The third point of Theorem 12.21 is obtained just by using the regularity result of 
v (point 2 of Theorem 12.2)1 . while in the second theorem we use in a crucial way the gravitational 
coupling Ap = p — 1. The techniques are thus different, and the one employed in the first theorem 
could probably be used for other types of coupling. 

Remark 4- The consistency with smooth solutions is detailed in Theorem 17.11 We try there to 
examine in what class our variational solution is the unique solution of the (E — P) system satisfying 
the boundary conditions. 

Remark 5. In many assertions we only state that "0 can be chosen in such a way that..."; this is 
because <f> is uniquely determined only in the dp a.e. sense. However, p is uniquely defined, and v is 
unique dp a.e.. 



3 Existence and uniqueness of a minimizer for the action 

This section is devoted to the proof of 

2d 

Proposition 3.1 Under the assumption that po an d Pt ore in L d + 2 , there exists a unique minimizer 
(p, J,Vp) in C([0,T];V(T d ) - w*) x (M{D)) d x L 2 (D) for the Problem\U& under the constraints 

f23 El El UW- 
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The proof of this result will be obtained in two ways: the first will be based and convexity / l.s.c. 
properties of the functional we are minimizing. However, despite the simplicity of this proof, we 
will introduce another way of obtaining the minimizing solution, by studying a dual problem similar 
to the the Monge-Kantorovitch problem, classically used in optimal transportation problems. This 
technique will then be a useful tool to obtain regularity results. Therefore we will show the following 

Proposition 3.2 Let K = K(pq,pt) be the infimum of Problem II., % and assume that K < +oo ; 
then 

K= sup \ <j>{T) dp T - 0(0) dp + / q- |Vg| 2 /2 dtdxl . 

dt<t> + q+ |V</»| 2 /2 < 

Remark. In the next paragraph, we are going to show that indeed K < +oo, assuming some 
integrability conditions on the initial and final densities. 

3.0.1 Existence of an admissible solution 

We are now going to prove that the infimum of Problem II. 31 is finite. This proof will use the results 
obtained in p] and j2l], concerning the time continuous formulation of the Monge-Kantorovitch 
transport problem. The optimal solution for the optimal transport problem with quadratic cost 
turns out to have a finite action for our functional /, provided the initial and final densities are 
sufficiently integrable. We will use the displacement convexity property obtained by McCann, in 
order to show that at intermediate times the density remains in H~ l {T d ). 

Lemma 3.3 Under the assumption that po, pr are in L^+^{T d ), there exists (p, J = pv,p) satisfying 

(EJ El El UM and such that Hp, j >p) iS fi mte - Moreover p G L°°([0, T\\ L^(T d )) and v G 
L°°(D,rfp). 

Proof. We use the following result that combines jH] and |2"Tj : 

Proposition 3.4 Let po and pt belong to V(T d ) H L k {T d ) for some 1 < k < oo. There exists a 
unique pair (p(t, x), J = pv(t, x)) that minimizes the action 

A(p,pv) = / \v(t,x)\ 2 dp(t,x) 

J[Q,T]xT d 

among all (p, J) that satisfy p G C([0, T]; V(T d ) - w*), J G {M{[0,T] x T d )) d and 

d tP + V ■ J = 0, 
p(t = 0) = po, 
p(t = T) = pr. 

A(p,pv) is finite, and for k > 1 — 1/ d, the function t — > (k — l)\\p(t, is convex for t G [0,T]. 
Hence for k>l, \\p(t, ■) ||* fc is bounded by max{||p ||^ fc , ||pr||^fc}- Finally v G L°° (D , dp) . 

The minimization problem stated here is one formulation of the classical optimal transportation 
problem with quadratic cost. It is described in more details in paragraph 18.2.21 
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Using classical elliptic regularity we have, for 1 < k < oo, \\D 2 A~ 1 [p(t) — 1] Hifcnrd) < || p(^) || z, fe (T d ) • 

The Gagliardo-Nirenberg inequality gives ||Vp||x,fc*(T<i) < Ck\\D 2 p\\ L k(jd) with k* = Therefore, if 
k = we have 

|| VA- 1 ^,-) - l]|| L2(Td) < C\\p(t,-)\\ Lk{Td) . 

It follows that (p, J = pv, p) is an admissible pair for which / is finite. This completes the proof of 
Lemma 13.31 

□ 



3.1 Proof of Proposition 13.11 

We present here a proof of existence and uniqueness of the minimizer of Problem 11.31 As we have 
already seen, the functional 

dp+- I |Vp| 2 dxdt 
2 Jd 

is convex. The second term is strictly convex in p (if we restrict to probability measures) and lower 
semi-continuous (l.s.c.) with respect to the weak-* convergence of p. To see that the first term is 
also l.s.c. with respect to the weak-* convergence of the measures (p, J), we refer to ^H] where the 
following result is shown: 

Theorem 3.5 Let (p n , / n )neN be a sequence such that for all n 6 N, p n is a Radon measure on D, 
and f n is dp n measurable. Assume that p n — p for the weak-* topology, with p a Radon measure. 
Assume that j D dp n \f n \ 2 is uniformly bounded. Then there exists a subsequence (still labeled by n) 
and a pair (p, f), with f dp measurable, such that 

1. for all if e {C°(D)) d , / dp n f n -(p^> / dpf ■ (p, 

Jd Jd 

2. j D d ^\f\ 2 <liminf J dp n \f n \ 2 . 

Considering a minimizing sequence (p n , J n ) for Problem II .31 we get also a sequence (p n , v n ) where 
v n is the density of J n with respect to p n , that belongs to L 2 (D, dp n ). Thus (p n , v n ) will be as (p n , f n ) 
in the Theorem above. 

We will show that the sequence p n is equicontinuous in C([0, T]; V(T d ) — w*). For this we use 
the mass conservation equation dtp n + V ■ J n = 0. For any time t, we can estimate the total mass of 
J n by Cauchy-Schwartz inequality: 

M(J n (t)) < ( f dp n (t,x)\v n (t,x)\ 2 
\Jj d 

hence J n G L 2 ([0,T}; (M(T d )) d ) C L 2 ([0, T}; H~ s (T d )) for s large enough. Therefore from the mass 
conservation equation, d t p n is uniformly bounded in L 2 ([0, T]; H~ s ' (T d )) form some s'. Of course, p n 
is also bounded in L°°([0, T}; M. (T d )). Using standard arguments of functional analysis (seepH]). we 
obtain that the sequence p n is equicontinuous in C([0, T}; V(T d ) — w*). 

The first point shows that a subsequence (p n , J n = p n v n ) will converge weakly to (p, J = pv). 
This implies that one can pass to the limit in the continuity equation d t p + V • J = 0. Moreover the 
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equicontinuity of the sequence p n yields that p(0) = po, p(T) = px- We thus obtain that the pair 
(p, J) satisfies (US IH Efl • 

The second point shows that I(p, J,p) is smaller than liminf I(p n , J n ,Pn), and therefore equal to 
K = inf /, since we consider a minimizing sequence. 

For the uniqueness part, since the term J | Vp| 2 is strictly convex with respect to p, two minimizing 
solutions must have the same density p. Once the density is set, the velocities must coincide dp a.e. 
again by convexity with respect to v of v — > p\v\ 2 . This proves Proposition 13.11 

□ 



3.2 Introduction of the dual problem 

We use here standard convex analysis arguments that can be found in [Oj and the proof is an adap- 
tation of the one found in j?j . The constraints (jlUl fTTl IT2"1 [T^j) can be formulated in the following 

weak way: 

(16) V0 eC°°(D), [ d t (f>{dp-dp)+V<f>-{dJ-dJ) = 0, 

Jd 

(17) Vg e C°°{D), [ (dp - dp)q = - I (Vp - Vp) • Vg dtdx, 

Jd Jd 

where (p, J, p) satisfies the constraints (|10j) to (|13j) . hence the triple of Lemma 13.0.11 works . Minimizing 
/ under the constraints (fT0|) to fT3|) is thus equivalent to find 



K = inf sup < / dpc + dJ ■ m — dt<p(dp — dp) — V0 • (dJ — dJ) 

P>J>P cf>,q,c,m [Jd 



J ^|Vp| 2 - Vg - (Vp- Vp) dtdx - J q(dp - dp) J 



with the supremum taken over all the continuous functions c, m with c : D — > R and m : D — > M. d 
satisfying c + \m\ 2 /2 < 0. 

Here C(D) is the space of continuous functions on D and C#(D) is defined by the additional constraint 
that the integral over T d vanishes for all t G [0,T] . On C(D) we have the usual duality bracket 

< f,g > denoted by / f dg with g 6 A^(-D) the set of bounded measures on D. The dual space of 
Jd 

C# is reduced to the set of bounded measures g on D whose total mass at any time is zero (i. e. for all 
g G C#, for all z G C°[0,T], f D z(t)dg = 0) and denoted by Ji4#(D). We introduce the functionals 
a and (3 defined on (c, m, r) G C(-D) x (C(D)) d x C#(D). It will be convenient to denote r = Ag, 
and this is possible since the mean value of r is zero. 



a{c,m,r) r> 



\ [ |VA _1 r| 2 dtdx = \ [ |Vg| 2 dtdx 
2 Jd 2 J D 



if c + \m\ 2 /2 < 0, 
a(c,m,r) = +oo otherwise; 



(3(c,m,r) = / dpc + dJ-m+ / pr dtdx 
Jd Jd 

if 3 G C^-D) such that c + 9 t + g = 0, m + V a 
(3{c,m,r) = +oo otherwise 
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with (p, J, p) as above. 



We compute a* and (3* the Legendre-Fenchel transform (see jH] for definition) of respectively a and (3. 
They are defined on (p, J,p) G M(D)x(M(D)) d xM#(D) the dual space of C(£>) x (C(D)) d xC # (D). 
We have for a 



a*(p,J,p) = sup < / dpc + dJ-m + / rp — | Vg| /2 dida; f , 

e+|m| 2 /2<0, r=Ag JD 



which is equal to 



<x*(p,J,P) = \J D 





j 


2 




L 




dp H 


-f 




p 


2Jd 



The term 



dp is defined through (JHJ); note that this can possibly be +oo. Then for (3 we have: 

)■ 

the supremum being restricted to all the c, m, r = Aq such that there exists <fi satisfying: 



f3*(p, J, p) = sup < / (dp — dp)c + (dJ — dJ) ■ m + / (p — p)r dtdx 

c,m,r {Jo JD J 



c + d t (f) + q = 0, 
m + V x = 0. 



Thus in terms of (f>, q we have 



sup < (dp — dp)(-d t (j) -q)- (dJ - dJ) ■ W x <p - / Vg ■ (Vp - Vp) rftdx 

Using the fact that (p, p) satisfy (JTDHHJ) we find that (3*(p, J,p) = if (p, J, p) satisfies ([TH | IT7 ]l and 
13* (p, J,p) = +oo otherwise. It follows that 



K= inf {a*(p,J,p)+(3*(p,J,p)}, 
p,J,p 

where we now compute the infimum over all (p, J, p). We have just relaxed the constraints 
ffTUl fm U2J O by adding the convex functional (3* which is +oo if they are not satisfied and 
if they are satisfied. 

3.2.1 The duality theorem 

The functions a, (3 are convex with values in ] — oo, +oo]. At point c = — l,m = 0,r = 0, 
a(— 1, 0, 0) = 0, a is continuous with respect to the norm of C(D) x (C(D)) d x C#(D), and 
(3(— 1,0, 0) = — J D P — —1 is finite. The conditions to apply Fenchel-Rockafellar duality Theorem 
(see ch. 1] ) are thus fulfilled and we obtain that 



inf {a*(p, J,p)+(3*(p, J,p)} 
sup{— a(— c, —m, —r) — /3(c, m, r)} 
if. 
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and the infimum is attained. Hence we have 



K= sup { / -dpc- dJ -m+ / Vp-Vq - \Vq\ 2 /2 dtdx} 

c,m,r=Aq JD JD 



c = -d t (j) - q, 
m = — V0, 
-c+ |m| 2 /2 < 0, 

which is equivalent to 

K = sup{ dp (d t <p + q) + dJ ■ V0 + I -| Vq\ 2 /2 + Vp ■ Vq dtdx}, 
<t>,q JD Jd 

the supremum being performed over all (<j), q) such that 

d t (f) + q+ |V0| 2 /2 < 0. 

To highlight the analogy with the Monge-Kantorovitch duality, and using (fTUllTTl [T^llT^j) . we rewrite 
the dual problem as follows: 

K= sup \ <j){T) dp T - 0(0) dp + / q - |Vg| 2 /2 dtdx) . 

The first integral is as in the Monge-Kantorovitch problem, and the second takes into account the 
non- linear effect of the gravitational coupling. This shows Proposition 13.21 

We choose p = p = 1 + Ap and J = pv to be any optimal solution (i.e. any minimizing solution). 
Note that necessarily J has a density v with respect to p and v G L 2 (D,dp). This justifies the 
notation J = pv. Then for all e > there exists (4> e ,p e ) G C l (D) with d t (p e + Pe + |V0 e | 2 /2 < such 
that 

K = \ [ dp\v\ 2 + l [ |Vp| 2 dtdx 
2 J d 2 J D 

(18) < / dp (dt<pe + Pe) + dJ ■ V0 e + / - 1 Vp e |7 2 + Vp • Vp e dtdx + 



e 2 . 



thus 



-/ dp\v -V0 e | 2 + - / |Vp- Vp e | 2 dtcfe 
2 J d 2 J D 

< [ dp(d t <p t +p € + \V^\ 2 /2) + e 2 , 
Jd 



and we obtain 



/ \dp\v- \7(j) e \ 2 + dp\d t <p € +p e + |V0 e | 2 /2| + / -|Vp- Vp e | 2 dtdx 
2 Jd 2 



<e 2 . 



It follows that, as e — > 0, 

• V0 e converges to v in L 2 (D, dp), 

• Vp e converges to Vp in L 2 (D, dtdx), 

• <9 t e + p e + I V0 e | 2 /2 converges to in L 1 (-D, (ip). 
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3.2.2 Uniqueness of the minimizer 

This property has already been proved directly, we just explain how to recover it from the dual 
formulation. Notice that the sequence (0 e ,p e ) does not depend on the optimal solution (p, v) we have 
chosen, thus if we have Vpi and Vp2 two optimal solutions, then Vp e converges to both Vpi and 
Vp2, and they are equal Lebesgue a.e.. It follows then that two optimal solutions have the same 
density p. Then since V0 e converges to both v\ and t>2 in L 2 (dp), v± and v 2 are equal dp a.e. which 
proves the uniqueness of the optimal solution. This ends the "dual" proof of Proposition 13.11 

□ 



4 Variations near optimality 
4.1 Outline 

The results of the following two sections require to establish the formula (J20|) . We present first 
the general approach of the proof, before entering into the rigorous details. We first state a simple 
property of convex functions near their infimum: Let \l/ be a convex function, and suppose it reaches 
its minimum value, say at x = Xq. Assume that D 2 ^> > XI. Then 



\x - x \ 2 < -r(^(x) 

A 



The formula (J20|) expresses this fact. Then bounding by the left hand side, we get that V\l/(xo) = 0. 
This yields the optimality equation as well as the conservation of energy. Keeping the left hand side, 
we can also estimate \x — Xq\ in terms of — This yields the formal H 1 (dp) estimate of 
Proposition 16. 11 In the case studied here, a formal calculation yields 



I(p2, h-,V2) ~ I{pi, JuPi) 



\Jo. 



P2 



+ |Vp 2 



pi 



+ |Vpi 



(19) 



+ 



■h 
Pi 
1 



{J2-Jl)-]~4-{p% 
z p l 



Pi) + Vpi • (Vp 2 - Vpi) 



:P2 



1 



|Vp 2 - Vpi 



J2 ^ Jl 

P2 Pl 

This equation is the formal analogue of formula ()20|). The second line contains the first order terms, 
and the third line is a positive function. Assume that (pi, Jx,px) is a critical point of /, the second 
line vanishes. First for J 2 = J\ + w where V • w = 0, we have p 2 = pi, this gives that J Vi • w — 
for all divergence free vector field w, thus v± is a gradient, which we already knew, thanks to the 
dual formulation. Then using that d t (p2 — Pi) + V • ( J2 — Ji) = 0, with p 2 = pi at times 0, T, using 
Ji/ Pi — V0i, and the Poisson equation Ap = p — 1, we get for the first order terms (second line) 



-9 t (f>i(p2 - Pi) 



{p2 ~ Pl) ~Pl{p2 ~ Pl) 







for all P2. This gives the Hamilton- Jacobi equation 



^0i + ^|V0i| 2 + Pi = 0. 

(Note that thanks to the dual formulation, we had an approximate "dp a.e." version of it.) Taking 
the gradient of this equation gives the momentum equation (JHJ) 



d t vi + vi ■ Vf 1 = -Vpi 



Vl = V0i. 
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This will be the Proposition 15.11 

Now at the critical point we keep only the third line of equation (JT9j) (the second line vanishes 
for all p 2 ), and we get 

-p 2 \v 2 ~ ^l| 2 + - |Vp 2 - Vpi| 2 = I{p 2 , J 2 ) ~ I(pt, Ji). 

If we perturb pi in p 2 = (x + 5(t,x))#pi for a small smooth perturbation 5(t,x), and perturb V\ in 
order to satisfy the conservation of mass (J2J), we will obtain an estimate of the form 

/]_ ]_ 
-Pi\v 1 {x + 5{t,x))-v 1 {x)\ 2 + - \V Pl {x + 5{t,x)) -V Pl {x)\ 2 

< C J \d t 5\ 2 . 

This estimate is a kind of if 1 (<ipi) estimate for v\, and also a if 2 estimate for pi, therefore a L 2 
estimate for pi. We will then adapt techniques from ^T] to obtain strong time regularity for p\. This 
will yield p 1 G C(]0,T[; L p (T d )) for p < 3/2. This will be the Proposition IO 

However, all these calculations are formal, since we have to integrate v 2 against the measure p±, 
and vice- versa, and we control V{ only in L 2 (D, dpi). Thanks to the dual formulation we will overcome 
this difficulty, by defining an extension of v in L 2 OC (]0, T[xT d ). 



4.2 Second variation formula 

We first introduce a perturbation of the optimal path. We must perturb the optimal pair (p, J) in 
such way that the conservation of mass (equation (flfi|)) is still satisfied. We proceed as follows: 

Let 5 and p be two small parameters and take r e]0, |-[. Let £(t) be a smooth function compactly 
supported for < t < T. We choose p small enough such that t — > t + pC(^) is a diffeomorphism 
from [0, T] to [0, T]. Let x — > w(x) be a smooth vector field and (s,x) e sw (x) the flow associated 
to w(x) defined by 

d s e sw {x) = w{e sw {x)) and e 0w (x) = x, 
we can thus define e l5 ^*- )u '(x). We introduce, as in [7], the following measures: 

p r >(t,x)=p(t + V ((t),x), v ri {t,x)=v(t + r 1 ((t),x){l + ri((t)). 

We check that the pair (p v , p^v 71 ) satisfies the continuity equation (|TE|) . We define also p n = p(t + r]((t),x). 
Then we define the measures (p v ' s , J v ' s ) so that for every / G C(D) and g G {C(D)) d we have 

/ f(t,x)dffi> fi (t,x)= [ f(t,eK® w (x))d/P(t,x), 

J D J D 

and 

/ g(t,x)-dJ v ' s (t,x) 
= [ g(t,e scit)w (x)) ■ [(dt + vifax) ■ V)e 5<:{t)w {x)}dp r '(t,x). 

Note that this can be rewritten We check that the pair (p^ 5 , J v > 5 ) satisfies also the 

continuity equation (|16|). Parameters p, 5 being fixed, we will use the following notation: 

v v '\t,x) = (d t + v v (t,x) ■ V x )e 5(:(t)w (x). 
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(Note that we have v v ' 5 (t,x) = v' n (t,x) + 0(5).) Therefore, J v ' s can also be rewritten e S %p n v r, ' s . 
Considering (0 e ,p e ) a smooth maximizing sequence for the dual problem, we have 



p%t,x)(d t cf> e + l ^l +Pe )(t,e 5 ^ w (x)) < 0, 



and using (JTSJ) we can write: 



\j D dp\v\ 2 + \Vp\ 2 



< e 2 + / dp(d t t + p e ) + dpv ■ V0 e + Vp • Vp e - | Vp e | 2 /2 
Jd 



- I dp^(d t <p e + ^^+p e ). 



Then using the mass conservation equation (fTp^) we have 

1 f dp\v\ 2 + |Vp| 2 < e 2 + / dpv (v^ s - V0 e (e^)-^|V0 e | 2 (e <5 ^; 



2 ./£) in V 2 



- dp v p e {e SCw ) + Vp • Vp e - |Vp e | 2 /2, 



so Ap = p — 1 yields 

^ / dp\v\ 2 +\Vp\ 2 < e 2 - \ [ dp^v"' 6 - VMe 5Cw )\ 2 + l [ dp^ s \ 2 
2 Jd 2 J D 2 J D 



I p e -p e {e 5 ^) 

D 



+ / D{e^ w ) : Vp v ® Vp e (e 5 ^) - |Vp e | 2 /2, 

where D(e l5< '"') denotes the spatial derivative of We obtain the complete formula: 

dp^\v^ 5 - Vcj) e {e^ w )\ 2 + \ ( \Vp v - D{e^ w )V Pe {e s < w )\ 2 
> 2 J D 

+ \ J D dP> V ' & \ 2 -\ J D d P\v\ 2 



< e 2 



+ ±JjD{e^)Vp t {e Kw )\ 2 -\J D \Vp e 

+ I f IW| 2 -~ f |Vp| 2 



2 



2 jd 2 j D 



(20) + I Pt- Pe (e s ^ 

D 



5 Optimality equation 

In this section we prove the following: 



17 



Proposition 5.1 The solution of Problem M.^M s a weak solution of the Euler-Poisson system (E — P) 
in the sense of Definition \2. l[ The energy of the system defined for a.e. t e [0, T] by 

E(t) = i / dp(t,x)\v(t,x)\ 2 - \Vp(t,x)\ 2 dx 
2 Jjd 

does not depend on time. 

Remark. The energy is a priori well defined in L 1 (0,T) since I(p,v,p) is finite. 

5.1 Derivation of the momentum equation ((31) 

In (I20J) . taking r\ — 0, bounding the L.H.S. from below by and letting e go to we get: 

0< \ JjPlv^-l f D d P \v? 



+ I^|D( e ^ )Vp (e^)| 2 -i^|Vpp 
+ f p-p(e s < w ). 



D 



Remark. Notice that this inequality (or a similar one up to second order terms in 5) could have 
been obtained directly by expressing that 

I(p 5 ,J S ,p 5 )>I(p,J,p), 

where p s = A~ 1 (p s — 1). Notice however that the expression of p s is not so straightforward. 

Expanding e s ^ w {x) = x + 5((t)w(x) + 0(5 2 ), we get 

\D(e s ^ w )Vp{e^ w )\ 2 

= \J W + KDw){e- & ^ w )Vp\ 2 J(e- 5 ^) + 0{5 2 ) 

with / the identity matrix of order d and J(e _<5 ^^- )u ') the jacobian determinant of the mapping 
x -> e~ s ^ w (x). Using that J( e - 5 ^ w ) = 1 - 5((t)V ■ w + 0(5 2 ), this is equal to 

\ [ (\Vp\ 2 + 25(Dw : Vp® Vp-5(\Vp\ 2 V- w) + 0(5 2 ). 
2 J D 



Then for v we have 



and 



~ J d *p\A*-\ J D d P\ v \ 2 

[ dpv (d t + v V)5(w + 0(5 2 
Jd 



- p ( e s ^) = - p- Vp-5(w + 0(5 2 ). 
d Jd Jd 
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This yields finally 



< 5 



dpv ■ (d t + v ■ V)(w + (Dw : Vp ® Vp 



D 



— -|Vp| 2 V • (w — Vp ■ (w dtdx 



+ 0(6 2 



Thus for every w smooth vector field on T d , £ G C£°(0, T), we have 

= dpv ■ (d t + v -V)((w) + (Dw :Vp®Vp 
Jd 

— -|Vp| 2 V ■ C w — ■ £w dtdx, 

hence we conclude that (p, v,p) is a weak solution of the Euler-Poisson system in the sense of Defini- 
tion 12.11 Note that in the proof of Proposition 13. 1[ we have already shown that 
p G C([0,T];V- w*). 

□ 



5.2 Conservation of energy 

Here we shall deduce conservation of energy by using a wrinkle in time to perturb the minimizing 
path; we take 5 = in (J20|) . minorize the LHS by and let e go to to obtain 

< I J dp{t + V ((t), x)(l + v((t)) 2 \v(t + V((t),x))\ 2 

+ \ I \Vp{t + r]({t),x)\ 2 dtdx 
2 Jd 

— - dp(t, x)\v(t,x)\ 2 - - / \Vp(t,x)\ 2 dtdx. 
2 Jd % Jd 

Changing variable in time t := t + rj((t), dt := dt(l +T]((t)) we get 

< \ I dp^x^v^x)] 2 ^) + \Vp(t,x)\ 2 (- - 1) dtdx. 

1 Jd 1 + i]Q{t) 

Taking the first order term in rj we get 

- / [dp(t, x)\v(t, x)\ 2 - \Vp(t, x) \ 2 dtdx] C(t) dtdx = 
2 Jd 

for any ( G C^°(0, T) which gives the conservation of energy, with 

E = - I dp(t, x) \v(t, x)\ 2 - \Vp(t, x)\ 2 dx. 

2 Jjd 

This ends the proof of Proposition 15.11 

□ 
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6 Regularity properties of the minimizer 



In this section we obtain several regularity properties of solutions of Problem 11.11 Those properties 
will follow from equation (f20|) . the rigorous analog of |T9|) . Here we use in a crucial way that our 
solution is a minimizer of the action of the Lagrangian and not only a critical point. For points 1 
and 2 we follow closely the method of Brenier in fJJ where similar results were obtained in the case 
of the Euler incompressible equation. For the fourth point we use the first two points and a method 
close to the one used by DiPerna and Lions in jllj . Similar results have also been obtained by 
different techniques in [Q] for finite dimensional hamiltonian systems, using properties of a special 
Hamilton- Jacobi equation related to the Hamiltonian flow. This section is thus devoted to the proof 
of the following: 

Proposition 6.1 The optimal solution (p, J = pv) of Problem M .S\ has the following regularity prop- 
erties: 

1. The density p belongs to L 2 oc (]0,T[; L 2 (T d )), 

2. the velocity v can be extended in all of T d to a function of L 2 OC (]0, T[; L 2 {T d )), in such a way 
that for all r in ]0, T/2), for all y in R d , 



3. the velocity potential (ft can be chosen in L^ c (]0, T[xT d ), 

4. finally p belongs to C(]0, T[; L p (T d )) for any p E [l,3/2[. 

6.1 Spatial regularity: proof of points 1 and 2 

We are going to deduce spatial regularity using the time dependent uniform translation 
e sc,w (x) = x + S((t)y to perturb the minimizing path. This corresponds to the case r\ = 0, £ = 1 in 
[r, T — t], w(x) = y fixed in inequality (J20|) . In this case v s (t, x) = v(t, x) + 5((t)y, and D(e sc - W ) = I, 
and inequality (|2l]j) becomes 





For the last line, we have used — dpv 

dt J T d 



0; indeed take tp = ((t)y in the momentum equation 



(JHJ), and use / Vp = 0. We have also 
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for a suitable choice of £. Hence 



I [ dp{t,x)\V(j) t {t,x + y) -v(t,x)\ 2 + \Vp e (t,x + y) -Vp(t,x)| 2 dtdx 

Jf d J : 

< e 2 + %\ 2 . 

T 

We let e go to and obtain 

f rT—T 

I I \Vp(t,x + y)-Vp(t,x)\ 2 dtdx<—\y\ 2 , 

Jf d Jt T 

thus D 2 p, and p = 1 + Ap are in L 2 OC (]0, T[; L 2 (T d )). In particular, p is absolutely continuous with 
respect to the Lebesgue measure of D. We will also obtain that V0 e is bounded in L 2 OC (]0, T[; L 2 (Y d )). 
Indeed, we get first from (J2TJ) that 

(21) / / p(t,x)\X7<p e (t,x + y)\ 2 dxdt<C{\ + -) 

Jj d Jt T 

for e < 1 . Integrating this over y G T d we get 

cT — T 



p(t, x + y) |V0 e (t, x)| 2 dxdydt 
[ [ \V^(t,x)\ 2 dxdt<C(l + -), 



y£l d Jx£T d Jt 
T-t 



thus we can, up to extraction of a subsequence, define a weak limit for this sequence, as e goes to 0, 
v = V0 G L 2 OC (]0, T[xT d ). However v will be uniquely defined only in the dp a.e. sense. Moreover 
for all y G T d , V0 e (- + y) converges to V0(- + in L 2 (dp) weak: indeed, from (}2*T]) . the sequence 
V0 e (" + y) is uniformly bounded in L 2 (dp), and for all <p G C£°(]0, T[xT d ), we have 

p«, x) V0«(x + ,) • „(t, x) *b / x)V0(x + „) • V (t, x) A*, 

since p G L 2 OC (]0, T[xT d ), and since V0 e converges weakly to v in L 2 OC (]0, T[xT d ). Using that the 
L 2 (dp) norm is l.s.c. with respect to the weak L 2 (dp) convergence, we will have 

P (t, x)| W(x + y)? dUx < U m inf / x,| V*(x + ,)|^x. 

Therefore -> V0(t,x + y) will be in L 2 oc (]0,T[; L 2 (cfe) n L 2 {dp{t))) for any y G T d . We finally 
obtain 

(22) f I p(t,x)\v(t,x +y) -v(t,x)\ 2 dtdx <C T \y\ 2 . 

Jf d Jt 

This proves the first two points of Proposition 16.11 
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6.2 L°° bound for the potential 0: proof of point 3 

We assume here that d < 3. We normalize p so that its mean value is for all time. Then 
since p G Lf oc (]0, T[xT d ) we have p G L 2 OC (]0, T[; H 2 (T d )) which is continuously embedded in 
L 2 oc (]0,T[; C5(T d )) and thus \\p(t, .)||l« e ^ OC (]0,T[). (We will see after that p is in L~ (]0, T[xT d ) 
and we will be able to remove this assumption on the dimension.) 
We take a regularization in t, x of <p,p: on ]r, T — r[, 

=T] e *(f), 
P e (t,x) =T] e *p, 

Ve(t,x) = -JTiVl(-,-), 

with 771 compactly supported in [—1, 1] x B(0, 1) , and < e < r/2. We first check that 

(23) ^ + i|V0f + p e <0. 
Considering cr — > </> e (cr, 7(0-)) with 7 G C 1 ([0, T]; T d ) we have 

^(0 e (a, 7 (^))) = ^(a, 7 (a))+7-V0 e ( ( x,7( ( T)) 

< ^(a, 7 (^)) + i|V0f (a n (a)) + ^(a)\ 2 

< -^(a, 7 (a)) + ^| 7 (a)| 2 

using (|2*3j) . and we obtain 

e (t + S ,x) < inf |^(t, 7 (t)) + ^ +S -^(a,7(a)) + ^|7r(cr)rfa}, 

with T the set of all continuous paths going from [t, t+s] to T d such that j(t+s) = x. Then restricting 
the infimum to paths of the form 7(0") = 7 (t) + ^f^ix — 7(f)) and noticing that 

lb e (*j OIU 00 e ^LcQOj^D (uniformly in e) implies that J* +s \p t (a,'j(a))\da < Cy/s, we obtain the 
following upper bound: 

(24) <p e {t + s,x)< inf {0 £ (t, z)} + C(~ + Vi). 

zeT d s 

A simple computation shows that 

(25) f dp(t 2 ,x)<j) e (t 2 ,x) - [ dp{t u x)4) e {t u x) 

Jj d JT d 

~* 1 2 I l\V<j)\ 2 dp+\Vp\ 2 dtdx 

Jti JT d ^ 

as e goes to 0, in particular it converges to a non-negative value, and will be greater than —1 for e 
small enough. It costs no generality to normalize <p e such that 

/ dp{T/2,x)<p e {T/2,x) = 0. 
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Then by choosing (£1,^2) = ( r /2,T/2) and (ti,t 2 ) = (T/2,T — r/2) in (j25j) . we get respectively that, 
for e small enough, 

(26) / dp(T/2,x)<f> e (r/2,x)<l, 



(27) / rfp(T-r/2,x)0 e (T-r/2,x) > -1. 
This implies 

inf{0 e (r/2,x)}<l, 

a; 

sup{0 e (T-r/2,a;)}>-l, 

and using (|2"4"j) we obtain 

(28) ||0 e |U~([r,T-r]xT«) < 

This leads to the conclusion, since e converges <it<ir a.e. to 0. 

Remark 1. This regularization of <j),p is needed to enforce the uniform bound on p e in H 2 . This 
might not have been satisfied by a smooth maximizing sequence of the dual problem. 

Remark 2: L 2 bound. Here we give also a bound when we do not assume d < 3. In this case we 
will have, from the Lf oc bound on p, p e G Lf oc (]0, T[; L k (T d )) for k = > 2. We then consider, 
for s,t g]0,T[, two probability measures p s and We consider the geodesic path between p s and 
p t , that we denote p s ',s f G [s,i], obtained from Proposition 13.41 (with proper time renormalization) . 
The corresponding velocity v(s'), s' G [s, t] satisfies d t p + V • (pv) = 0. Combined with (J23|) . this will 
yield 

<j) e (t)dp(t) < [ (j) t (s)dp(s) + / / (i|y(s / )| 2 -p e (s , ))rf/i(s')c/s'. 



Using the properties of the geodesic path p expressed in Proposition 13 A\ we are able to obtain that 



(j> e {t)dp:{t) - / ^(s)dp(s) <C(r,\t- s\)su V {\\p(s)\\ L i,\\p(t)\\ L2 }, 

with C(r, |i — s\) uniformly bounded for r > r , r < s, t < T — r, t — s > 5 > 0. (We have used in 
particular that |u| < C(d)/\t — s\.) In (J2TH I27|) above, we choose r such that p(r/2), p(T — r/2) both 
belong to L 2 with norm bounded by some constant C(r). Since p G L 2 OC (]0, T[xT d ), this is true for 
almost every r > 0. 

From this, taking successively p(s) = p(r/2),/i(t) = p(T — r/2), we obtain that for all probability 
measure p, for all t G [t,T — r], | f e (t)ci/i| < C(r)||yu|| L 2. This then yields 

(29) ll0llL°c([r,T-r];L W ) < C(r). 

6.3 Strong time continuity of p: proof of point 4 

First let us notice that in inequality (j2T)j) . if we do not set rj = 0, we obtain the mixed derivative 
estimate: 

-T-t r 

|2 



(30) / / P (t,x)\v(t + V ,x + y)-v(t,x)\' 

+ |Vp e (t + ^, x + y) — Vp(t, x)\ 2 dtdx 

< c T (\y\ 2 + \v\ 2 )- 

By taking y = we get that <%Vp G L 2 oc (]0,T[; L 2 (T d )). This implies 
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Lemma 6.2 The gravitational field Vp belongs to C£ c (]0, T[; L 2 (T d )). 

Now we prove the last part of Proposition 16.11 We obtain that the density is strongly continuous 
with respect to time by showing some renormalization property, that gives the continuity of some LP 
norm, with p > 1. Note that in the non-interacting case, since the functions of time (k — 1) \\p(t, ■) \\ k Lk 
are convex, they are continuous in ]0, T[ provided they are finite at t = or t = T. So we will prove: 

Lemma 6.3 Let a G [1, 3/2[ and G a (t) = / p a (t,x) dx, then G G C(]0,T[). 

We postpone the proof of the lemma after the proof of the last point of Proposition 16.11 

Proof of the last point of Proposition 16.11 First we check the weak time continuity of p: from 
the conservation of energy and from Lemma [6.21 J Jd p\v\ 2 is uniformly bounded on [r, T — r]. Thus 
P\ v \ = a/Pa/ p\v\ 2 G L°°([t, T — t];LP) for some p > 1 thanks to Lemma 16.31 It follows that from 
equation (J2J), d t p is bounded in L°°([t, T—t], LL~ s ) for some s. Using classical arguments of functional 
analysis (see [201) we can deduce that, for some p > 1, 

P eCQ0,T[,LP-w). 

Then Lemma 16.31 implies that p G C(]0,T[; LP) for any p G [l,3/2[: indeed it is a classical fact 
that when a sequence (« n ) n eN converges weakly in L p , 1 < p < oo, to some u, if ||w n ||i/p converges to 
|| u || lp the sequence converges strongly. The last point of Proposition 16. II is thus proved. 

□ 

Proof of Lemma 16.31 Let us prove the renormalization property when the density and the velocity 
field are smooth: we use the identity 

d t [pF(p)} + V ■ [pF(p)v] = - P 2 F'{p)V ■ v. 

Integrating over T d we get that J Jd pF(p) is continuous with respect to time as long as p 2 F\p)V ■ v is 
in Ll oc ([0, T] x T d ). We will see that this is true for F(p) = p@, /3 < | from the regularity property ([22)1 . 

I 1 2 

We introduce r](x) = Cexp( , ) with C such that L d r](x)dx = 1. Then as usual 

yl+|a;| 2 JK 

fj e (x) = -^T)(-), and p e (x) = p * r] t (x) = f Rd p(x — y)r] e (y)dy, p being naturally extended to a Z d 
periodic function on all of M. d . With (p,v) as before, we consider the pair (p e ,v e ) defined by 

Pe = Ve* P, 

v e = i] e * (pv)/p e . 

We check that the pair (p e ,v e ) still satisfies the mass conservation equation (|2j). Then we have the 
crucial property: 

Lemma 6.4 For r g]0, T/2[, 

Pe\V-v e \ 2 <C(r). 
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This lemma means that the spatial regularity property (}22|) is conserved through regularization. 
Before proving this lemma, we conclude the proof of Lemma RHIl since p t > 0, for all F G C' 1 (]0, +00 [) 
we have 

j t J p.F(p,) = -J p>,F'(p,)V-v„ 
|/rf=(l- a )/rfV.«, 
Then using Lemma ffi. 41 and the fact that p € G L 2 OC (]0, T[xT d ) we get that 



[ T T [ P 3 J 2 \V-v e \<C( 



r), 



and also that for any a e]l/2, 3/2[, the sequence pf\ V ■ u e | is equiintegrable on [r, T — r] x T d . Thus 
the sequence t — > f Jd pf(t,x)dx is equicontinuous and the limit t — ► L d p a (t,x)dx is continuous. 

□ 

Proof of Lemma 16. 41 

We have 



Vp e , . 1 , . _ 

— 5- • {pv) * Ve H (/w) * Vl] e 

Pe Pe 



Pi 



p(x - y)(v(x - y) - v(x))r) e (y) 



Vpe 

Pe 



^7 1 f 

H 1 • v H / p(x - y)(u(x - y) - f (x)) • Vr] e (y). 

Pe Pe JM. d 

We use the special shape of the regularization kernel: there exists C > such that for all x G M. d , 
\Vr)(x)\ < Crj(x). This implies 



< 



c 



We also have the usual bounds 



We define 



Then 



A{x) 



A (x) < C 



\y\\Vve{y)\dy<C, 

\yW(y)dy < Ce 2 . 



p(x - y){v(x -y)- v(x))r] e (y) dy. 



p(x - y)ve(y) \v{x -y)- v(x)\ 



Pe{X) 



dy 



and by Jensen's inequality this is less than 



r , , P( x -y)Ve(y)\v(x-y)-v(x)\ 2 

C I T-r dy, 

p e [x) e 2 
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hence we obtain 



v(x + y) — v(x)\ 



2 



p e (x)A (x)dx < / p(x)rj e (y) dydx. 

J d Jl d JR d e 



B(x) = — I p(x- y)(v(x - y) - v(x))Vr) e (y) dy, 



For the next term 



Pc 

we proceed by the same method. We now have the bound 

oT-T 



I I p e \ V ■ v e \ 2 (x) dxdt 

Jt Jf d 

T-r r r 1 



< C / / p(t,x) — \v(t,x + y)—v(t,x)\ r) e (y) dxdydt 

Jt Jj d JR d e 

= c I II P(^ x )ri9\ v ^ x + y) ~ v (^ x )\ 2 Ve(y)^- dxdydt 

Jt Jj d Ju d \y\ e 

< C(r) 

using the spatial regularity property ([22)1 . The proof of Lemma f6. 41 is complete. 



□ 



7 Consistency with smooth solutions of the Euler-Poisson 

system 

Here we show that the solution of the minimization problem coincides with a smooth potential 
solution of the Euler-Poisson system satisfying the boundary conditions, when the latter exists. We 
will try to make precise the smoothness required to reach this conclusion. 

Theorem 7.1 Let (V<f),p,p) be the solution of Problem \1.3l and let (Vip,r,q) be such that 

^ £ W 1,co (D), r E L 2 ([0,T];H- 1 (T d ))nL co ([0,T};L 1 (T d )), q G L°°(D). 
Suppose that (Vip,r,q) is a solution to 

(31) d t iP+^\ViP\ 2 + q<0, 

(32) r (^ + i|V^| 2 + g) = 0, 

(33) d t r + V • (rVV>) = m v '» 

(34) r\ t =o = pa, r\ t=1 = pt, 

(35) Aq = r - 1 in V. 

Then p = r and V0 = Vip holds dp a. e. . 

Proof. Since (p, V0) satisfies the continuity equation (|2j) we have 

pV0-W>+ / pd t i)= [ p T ip{T)- I Po ip(0). 

D Jd Jj d Jj d 
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Combining with (j31j) this yields 



f pV0-VV>- / P {\\V^\ 2 + q)> [ PtHT)~ 
JD JD 1 JT d 

Using the Poisson equation (jlj) we have 



Po^(0). 



1 y p(-|V0 - VVI 2 + | V0| 2 ) + J Vp-Vq-q 



> / p T 7/>(T) - / p ^(0) 

JT d JT d 

= / rV^ ■ + / rdtip 

JD JD 

= - J D ^ r+ \ y^iwf 

= J D \Vq\ 2 -q+ l -Jj\V^\\ 

where we have used (I33[ l34|) in the third line, (J32|) in the fourth line, and ()35)1 in the fifth line. We 
finally get 



/ p|V0-V^| 2 + / \Vp-Vqf 

J D JD 



< / (p|V</>| 2 + |Vp| 2 ) - / (r|VVr + |Vg| 2 ). 

Since (p, V0,p) is solution of the minimization problem, the RHS is non positive, and we obtain the 
expected result. 

Remark 1. This is true in particular if (ip, q) is a C 2 (D) x C X {D) solution of d t ip + || Vip\ 2 + g = 
and thus shows the consistency with smooth solutions of the (E — P) system. 

Remark 2. From the results of Theorem 12.31 depending on p and pT the assumptions on ip,r,q 
can be weakened. For instance if po an d Pt are in L°° then p, V0, Vp are in L°°([0, T] x T d ) and one 
only needs (1 + r)|VV ; | 2 , (1 + r^c^l, l^?| 2 t° be integrable to perform our computation. (Note that 
these assumptions imply that if} G C([0, T]; L 1 (T d )) and thus J Td pT^iT) — J Jd poV'(O) is well defined.) 
This ends the proof of the Theorem 17.11 

□ 



8 Path regularity for a time discretized interaction: Proof 
of Theorem E31 

In this section we prove several additional regularity properties for the variational solution. The 
problems that we will treat are closely related to viscosity solutions of Hamilton- Jacobi equation. In 
the remainder we will denote by HJ1, HJ2 the following operators: 

HJ1 (0)=^+i|V0| 2 , 
HJ2 (0, p ) = ^0+^|V0| 2 +p, 

and we will make precise the sense in which they must be understood. See J2| for references about 
Hamilton- Jacobi equations. 
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8.1 Formal bounds 



Solutions of our variational problem satisfy 

(36) d t <j>+ ^|V0| 2 + p = dp a.e., 

(37) ^+^|V0| 2 + p<O, 

(38) d tP + V ■ (pV0) = 0, 

(39) p = 1 + Ap. 

If some C 2 (T d ) function Z satisfies Z(xq) = and Z < in T d then D 2 (Z)(xq) < in the sense of 
matrices and in particular this implies that AZ(x ) < 0. Using (JBlfl l3T|) and applying this result to 
Z = d t <f) + \ \ V0| 2 + p one formally obtains 

d t A<f) + (V0 • V)A0 + \dij<P\ 2 + Ap < dp a.e., 

which combined with the inequality |A0| 2 < d^2 i - \dij4>\ 2 and with (|3T?j) gives the inequality 

(40) ^<l-p-||A*|", 

where the operator — • = d t ■ + (9,0 9j- denotes the convective derivative along the flow generated 



=i 



by the velocity field V0(t,x). We first deduce from this the following upper bound for A0 dp a.e.: 

A0<C(d)(l + ~), 
it is obtained by looking at the behavior of the differential inequality 

for large /. This "Oleinik-type" estimate is well known for viscosity solutions of HJ2=0 provided 
that p (as is the case here) satisfies Ap > C. However we don't know a priori that our solution is a 
viscosity solution, and moreover this bound is true in the sense of distributions. Here a complication 
is added by the fact that the solution (p,p satisfies HJ2(0,p) = only dp a.e., and therefore a priori 
not in the viscosity sense. 

Now notice that our variational solution exists on t e [0, T] and thus (-0, <i){t) = (— <p,p)(T — t) is 
also a solution to equations (|3H|) to (|3~9"j). Therefore we can obtain by the same way that 

(41) Am<C(d)(l + - t ) 1 

(42) A 0(T-t)>-C(rf)(l + i). 
This gives the following uniform bound 

l|A0|| L o°([r,T-T]xTd, dp) < C(d)(l + -). 
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Here the constant is universal, and it is only supposed that the solution exists from t — to t — T. 
This surprising (in the sense that it is not true for viscosity solutions) result comes from the fact 
that the transformation (<p,p) —>■ (—<f),p)(T — t) does not necessarily transform a viscosity solution 
of (J36|) into another viscosity solution but transforms a variational solution into another variational 
solution. Actually it will be proved that one can choose the variational solution to be a viscosity 
solution in one time direction, but it may only be a subsolution when reversing the time. However, 
it will coincide dp a.e. with the viscosity going from T to 0. This means that we have a reversibility 
property dp a.e.. 

Moreover A0 is the divergence of the velocity field, and we have A<fi = — J^logp. From ()40|) we 
have the following control on the second time derivative of p along the flow: 

(43) ^ logp > p -l + l|l, ogp | 2 . 

Following the path of a "particle", the differential inequality satisfied by = log p is: 

6 > exp6 - 1 + ^-lei 2 . 

d 

We look for solutions of this equation that do not become infinite in ]0, T[. This condition implies 
that 0(t) < C(r) for r < t < T — r, independently of the initial and final values of 0. We have thus 
an interior unconditional bound for the L°° norm of p, namely that 

IMU°°([t,t-t]xT<*) < C(r), V0 < t < T/2. 

The above differential inequality will also yield that some functionals of p are convex along the 
displacement induced by our variational problem: indeed a formal computation gives the following: 

d 2 f , 

— / plogp > 0, 



dt 2 
dt 2 



p\ k > for every k > 1. 



d 



Remark. This displacement convexity property is analog to the one found in [21] which was true 
for k > 1 — 1/d (when multiplied by k — 1). In our case this is only valid down to k > 1 due to the 
gravitational term. However, note that this displacement does not induce a distance: indeed take 
Po = Pt ^ 1 and check that the cost of the transportation of po on pt following the Euler-Poisson 
flow is not 0. 

In the next section we give a rigorous sense to the computations made above in order to obtain 
the Theorem 12.31 



8.2 Rigorous proof of Theorem 12.31 

Outline of the proof The rigorous justification will be achieved in several steps. We will first 
introduce a time discretization of the problem, where the potential energy term contributes only at 
the time ti = iT/N, i = 1..N — 1. Between two time steps, the problem will be an optimal transport 
problem as in jHj and j2l]. Then at each time step, the gravitational effect will be taken into 
account, and the velocities will have a "jump". From a Lagrangian point of view, the velocity of each 
particle will therefore be a piecewise constant function with respect to time. This discrete problem 
will accept also a variational formulation. Then letting the time step go to 0, we will eventually 
recover the time continuous problem. One of the advantages of this discrete formulation is that 
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between two time steps, the velocity potential will be expressed with convex functions, which offer 
several regularity properties, and that will allow us to deal rigorously with quantities such as the 
second derivatives of (p. Moreover, we will show that we can choose the optimal <fi to be the viscosity 
solution of the Hamilton- Jacobi equation HJ2(0,p) = 0. Then the Hopf-Lax representation formula 
will be a precious tool. 

This section is organized as follows: we first introduce the discrete problem, and recall some facts 
about optimal transportation^. e. the non-interacting transport). Then we show that solutions of 
this discrete problem converge indeed to the time-continuous one when the time step goes to 0. The 
rest is devoted to the proof of the regularity properties. For this we show that we can choose a special 
solution of the discrete problem that will be a viscosity solution. For this solution, we are able to 
perform rigorously our computations and conclude. 

8.2.1 Construction of a sequence of approximate solutions 

We introduce the discrete times ti —Ti/N, i — 1..N — 1, and consider the functional 

N—l 

I N (p,v,p) = ^ J^dp(t,x)\v(t,x)\ 2 + J^lVp^x^ 2 dx. 

We are now interested in solving the following variational problem: 
Problem 8.1 Minimize 

JV-l 



I N(p,J,p) = ^r T y2f \Vp(ti,x)\ 2 dx 
2/v ~[ Jr d 



+ sup < / c(t,x)dp(t,x) + m(t,x) ■ dJ(t,x) 

c,meC°(D)x(C°(D)) d ^ D 
c+|m| 2 /2<0 

among all (p,J,p) that satisfy p G C([0, T]; V(T d ) -w*), J G (M(D)) d , Vp(U) G L 2 (T d ) for all 
1 < i < N - 1, and 

d t p + V-J = 0, 
Ap = p - 1, 
p(t = 0) = po, 
p(t = T) = p T . 



We denote the value of this infimum. 

The interest of studying this problem is both to make rigorous the arguments of the previous section 
and to give a possible numerical discretization of the Problem 11.11 It will also let appear some 
interesting links between optimal transportation, viscosity solutions of Hamilton- Jacobi equations, 
and transport equations. 

8.2.2 Basic facts on optimal transportation 

We first recall the definition of the push-forward of a measure by a mapping: 
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Definition 8.2 Let po one? p\ be two probability measures on T d and let X be a dpo measurable 
mapping from T d into itself. We say that pi is the push-forward of p by X , that we denote by 
Pi = X#p , if the following holds: 



V/ G C\T d ),J f(X(x))dp (x) = J f(y)d Pl (y). 

We denote here Di = [tj, ti + \] x T d ; the effect of the time discretization is that between two times ti 
once the measures p(£j) and p(tj + i) have been chosen, the problem becomes the following: 



Problem 8.3 Minimize 



C(p,J) = sup {- / <I/>c -h (11 ■ in } 

c,m&C°(D l )x(C°(D,)) d 



D, 



c+|m| 2 /2<0 

among all (p, J) that satisfy p G C([tj, t i+ i\; V(T d ) — w*), J G (M(Di)) d , and 

d tP + V ■ J = 0, 

p{t = U) = pi, 
p(t = U+i) = Pi+i- 
The infimum is denoted C_(pi, Pi+i, \U — t i+ i\). 

Remark. Performing a dilatation in the time variable we see that C_(pi, Pi+i, t) = \C_(pi, p%+x, 1). The 
Wasserstein distance (of order 2) between p, and Pi+\, denoted ^(pj, Pi+i), is given by 

[W 2 (pi, pi+i)} 2 = C(pi, p i+ i, 1). 

This problem has been solved in j^j, j^j where it is shown that there exists a unique solution 
(p, J = pv) (v is only unique dp a.e.) that satisfies: 

v(t = U,x) = (V<f(x) — x) dpi a.e., 

ti+i — U 

d t {pv) + V- {pv® v) = 0, 

det D 2 (p(x)p i+1 (V(p(x)) = pi(x), 

with if a convex function. The second equation means that the particle move with constant speed. 
The third equation is the Monge- Ampere equation that is satisfied in the following weak sense: 

*i0 / , npd\ 



V/ G C\T d ),J f(V<p{x))dfH(x) = j f(y)dp i+1 (y). 

This means that V<p pushes pi forward to p%+i- The Wasserstein distance between two probability 
measures can be defined equivalently in the following ways: 



Wj(po,Pi) = inf / dp\v\ 2 /2 

P' v </[0,l]xT d 



sup / dp (x)(\x\ 2 /2-$(x)) + d Pl (y)(\y\ 2 /2-y(y)) 

$(x)+^>(y)>x-y Jf d 

= inf / ~\x — m(x)\ 2 dp 

m # po=pi Jj d 2 



sup / dpi(x)(j)(l, x) — dpo(x)4>(0, x) 

d t cf>+\X7cf>\ 2 /2<0. 
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where the first infimum is taken over all the pairs (p, v) satisfying 

dtp + v ■ H = o, 

p\t=o = Po,p\t=i = pi- 

Note that the last formulation is strongly reminiscent of the formulation of Proposition 13.21 Under 
some assumptions of absolute continuity of po, p\ with respect to the Lebesgue measure (see j2Hj 
for a complete reference), those three problems have a unique solution. For the first it has already 
been mentioned above. Then the optimal m is equal to V<p, the optimal pair ($, is equal (up 
to a constant) to (<p, <p*) with <p* the Legendre transform of <p (see the definition (|48|) below), and 
0(0) = $ — |x| 2 /2, (/>(1) = \x\ 2 /2 — \E'. For additional references about the Wasserstein distance the 
reader can also refer to j2l] and |22j . 



8.2.3 Existence of a minimizer for the approximate problem 

Following exactly the same method as in the first problem we can show the existence of a unique 
minimizer to the Problem 18. II In this way we obtain the following: 

Proposition 8.4 There exists a unique p^ and a dp^ a.e. unique = V0at solution of Problem 
\8.1l Moreover it satisfies: 

1. There exists C such that for any < r < T/2, — \\p N (ti, ■)\\L^<T d ) < — ■ 

T<U<T-T 



2. The solution (pn,vn = V0jv) is a weak solution of 



rp N-1 



d t (pNV N ) + V • (p N V N <8> V N ) = -PNJj ^2 S t=tiVPN(ti), 

i=l 

d t pN + V • {p N v N ) = 0, 
Apjv = Pn - 1. 



where the product pn^Vn is defined as in \2.1\ The pair {(Pm.Vn) satisfies 

N-1 

2' ' 1 ■ N 



1 T Af_1 

dt(f>N + ^|V0at| 2 + — <Wiv < 0, 



i=l 



and also satisfies for any U,l < i < N — 1 

T 

(f) N (tf,x) -4> N (ti ,x) < —=p N (ti,x) dx a.e., 

T 

4> N (tf,x) - ^> N (ti ,x) = -—Pn{U,x) dp N (U) a.e.. 

Proof. The proof is the same as the time continuous version therefore we will only sketch it briefly. 

Concerning the existence of an admissible solutions for Problem 18 . II note that we don't need that 
either po or px is in any LP since two probability measures on T d are always at finite Wasserstein 
distance and thus one can exhibit an admissible solution by transporting po on p(T/N) = 1 between 
t = and t = T/N , then letting p(tj) = 1 for i < N — 1 and transporting p(^^-T) on p?. Note also 
that the solution of the continuous problem is admissible for the discrete problem. 

Then having chosen an admissible solution (p, J = pv,p), the problem Jjy admits a dual problem: 
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Proposition 8.5 Let be the infimum of Problem \8.1\ Let 

D N (?P,q) = / dpdtip + dJ -Vi) 

J[0,T]xT d 
N-l 



+ TtYI I dp(t i )q(t i ) + (Vp(t i )-Vq(t i )-\Vq(t i )\ 2 /2)dx. 



Then 



K n = inf In = supDpf, 

P> J ,P f,q 



where the supremum is taken over all pairs (ip,q) € C X (D) such that 

N-l 

N 



^^|W| 2 + ^ J> = ^)<0, 

i=l 

moreover the infimum is attained. 

Remark. The dual functional can be rewritten as 

N-l 

I ! 

DN(i>,q) 



f dp{T)^{T) - dp(owo) + ^ \ q{u)- l -\v q {u)\ 2 dx. 



The proof of this proposition uses the Fenchel-Rockafellar duality Theorem as the time continuous 
one. The optimal velocity will here also be potential. Then taking (p,v,p) = (pn,^(Pn,Pn) the 
optimal solution, for any e > we find if) € , q e such that D N (ip t , q e ) > K N — e 2 and we obtain: 

i r . ir w 

(44) 



\ [ dp N \\7<P N -ViJ e \ 2 + ~y2 [ \Vp N {U)-Vq e {t t )\ 2 dx 

r i t N ~ l 

+ / dp N \d t ij e + -\Vij e \ 2 + -J2$t=tMU)\ <e 2 . 
J D i=i 

Then perturbing (p N , v^) as we did in the time continuous problem, we find the optimality equation 
and the regularity properties. 

Now we prove the last 2 points of Proposition 18 A\ that link <f)(t~),<f)(tf),p(ti). Before this we 
point out that the L^ c (]0, T[; L 2 (T d )) bound for established in Remark 2, paragraph 16 .2\ is still 
valid for <pN with minor adaptations. 

Then, in every ]tj, t i+ i[, d t 4>N < ~ |V0jv| 2 /2 is a negative measure, in the sense of V'Qti, t i+ i[xT d ). 
In particular, for a.e. x, (f>N is decreasing with respect to time. It follows from the monotone 
convergence theorem that t — > <pN{t, •) has left and right (strong) limits in L 2 (T d ) at every t e]0,T[. 
This gives sense to (f)(t~) , <f>(tf) . 

Since 4>n is a subsolution to d t (p + | V0| 2 /2 + 5tj){U)T / N = on U+i[, we have immediately 
that <f>{tf) - <j){tT) < -p(U)T/N. 

Then we take the maximizing sequence ip e , and if a smooth function on T d , for which we have 

rti+S r rU+5 f 1 T f 

/ / pdebw = / / -P-|V^|V-T7/ p{U)<P + 0(e 2 ) 

[Pn{U + 5)ip e (ti + S) - p N {U - 5)ip e (ti - 5)]<p 

T d 

k+5 [■ 

p N v N ■ V(ip e (p). 

i 
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U-5 JT 



For the first equality we have used equation (jB|) and we have used the mass conservation equation 
for the second one. 

The functions ip t (t,x) are decreasing with respect to time in every and converging in 

L 2 oc (D) — w to 07V; therefore if) e (t) — > 4>N{t) in L 2 — w for almost every t e]tj, U + i[ (actually at every 
t where where 0/v is weakly continuous in time, thus almost everywhere, since 0/v is decreasing). 

Hence i/j e (t±5) — > 07v(t±<5) weakly for almost every S. As 5 — > 0, p(ti±5) converges weakly to p(ti) 
in L 2 : indeed, note that pN{ti) £ L 2 for 1 < 1 < N — 1, and between two t iy the problem coincides 
with the optimal transportation; then from Proposition I3.4[ we get p^ G C(]0, T[; L 2 (T d ) — w)). 

Then, since 0/v is decreasing with respect to time and bounded in L 2 , we get that 4>n{U ± $) 
converges strongly to 0/v(*i ) in L 2 as 5 goes to 0. The other integrals go to when 5 goes to 0, 
except the one involving p(ti), and we obtain 

T 

<M*+) - <M^ ) = -j^p N {ti) dp(ti)a.e.. 

8.2.4 Regularity properties of solutions of the discretized Problem 18.11 

Here we state the main result of this section, from which Theorem 12.31 will be deduced. 
Proposition 8.6 Let (t,x) — * (4>N,J n = Pn^(Pn,Pn) be solution of Problem \8.1\ 

1. There exists C depending only on T and on the dimension such that for all t in ]0,T[, 0/v(i) 
is dpjsr(t) a.e. twice differentiable and satisfies 

-C(l + — !— ) < A<p N (t, .) < C(l + h dp N (t) a.e.. 



2. The density pN is bounded in L^ c (]0, T[xT d ) uniformly with respect to N and belongs to 
f] k>1 C(]0,T[;L k (T d )). 

3. There exists C such that for any 1 < k < oo 

+ \)<j t log OIUh*)) < c(i + j^). 

I Tke function, [ lprl{t , x)] - dx , ( PnM , X ) dx are un.fonnly UpscMt, m tk respect to 

Jf d Jl d 

time in every interval [r, T — r] for t e]0, T/2], and converge as N oo to convex functions 
on [0,T]. 

5. The potential 4>n can be chosen such that for every < r < T/2, 

||V0/v|| L ^([ TiT _ T ] xT d) < C(r). 

6. One can also choose 4>n to be the viscosity solution of d t <pN + ||V0/v| 2 + St=tiPN = 
in the sense of < f5^) . 



7. All these results and bounds do not depend on p neither on pr, and are uniform with respect 
to N. 

The proof of this proposition is postponed to the end of the paper. First we use it to show the 
convergence of solutions of Problem 18. II toward the solution of Problem ll.il 
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8.2.5 Convergence of the solutions of the discretized Problem 18.11 to the solution of 
the continuous Problem 11.31 

Proposition 8.7 Let (p N , vn,Pn) be as above, and (p,v,p) be solution of the minimization Problem 

. . 2d 

with same initial and final densities in L d + 2 , then 



lim / / (p\vn — v \ 2 + | X?pn — Vp| 2 ) (t, x) dxdt =0, 
moreover pn^4>n converges strongly in Lj oc (]0,T[xT d ) to pV0. 

Then the last two propositions combined will yield the Theorem 12.31 when passing to the limit. 
Proof of Proposition 18.71 

Here we prove a slightly weaker version of Proposition 18.71 that allows us to pass to the limit in 
Proposition 18.61 and obtain Theorem 12.31 Then we can get the full Proposition 18.71 We first choose 
t e]0, T/2[ and k = k(r, N) such that < t k -i < r < t k < .. < t N - k <T - t < tu-k+i- We set 



F T M = \ ! f dp\v\ 2 + \ [ [ \Vp\ 2 dxdt, 

F T N {p,v) = \f [ dp\v\ 2 + ^-J2f \Vp(ti,x)\ 2 dx, 
1 Jo Jr d Z1 * ~1 J^ d 



understood that p satisfies Ap = p — 1. We need to introduce those truncated functionals since we 
don't know a priori that the potential energy remains bounded near the boundary of the time interval 
and thus the convergence of the Riemann sum to the integral is not clear. We shall see after having 

j——. 2d 

proved the Theorem 12.31 that this is the case when p and pt are in L d + 2 . 

It follows from Proposition 18.41 that there exists a unique minimizer (pn,vn — '^'Pn-iPn) for the 
functional under the constraints of Problem 18. II It can also be checked later in the proof that the 
regularity results of Proposition 18.61 remain uniformly valid for r < To. We consider (p,v) solution of 
Problem 11.31 From Lemma l6l2| ||Vp(i, Oil x, 2 rir rf ) ^ s con ti nuous in ]0, T[ thus 

Vr > 0, F T N (p, v) -> F T {p, v) as iV -> oo. 

Then F T (p,v) — F(p,v) — > when r goes to 0, and there exists a sequence (t n ) N£N * with I 
such that 

F T N »{p,v)^F{p,v). 

In the remainder of this proof we set Fn := Fjf (therefore r, k will both depend on N) and (pn, vn) 
the minimizer of Fj^ . Thus 

Fn(p,v) -> F(p,v). 

The idea of the proof is to show that F(pn,vn) — > F(p,v), and to use the coercivity of F to 
deduce that (pn, v n) ~^ (p, v )- The first point will need some regularity property for (pn, v n)- 
We have first F N (p N ,v N ) < F N (p,v) and lim sup F N (p N , v N ) < F(p,v). 
We claim that we also have lim |Fjv(pjv, vn) — F(p^, vn) \ = 0, this will imply 

lim F N (p N , v N ) = lxmF(p N ,v N ) = F(p,v) = lim F N (p,v). 
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2d . . 2d 

Indeed, since po and pr are in L d + 2 , point 3 in Proposition 18.61 implies that pn is in L°°([0,T]; L d + 2 ). 
It follows that the sequence Vp^ is uniformly bounded in L°°([0, T]; L 2 ) since Ap^r = Pn ~ 1 and 
using Gagliardo-Nirenberg inequality (cf. proof of Lemma f3.0.1|) . 

Then, we have d t f | Vpat| 2 /2 = — f Pn^n • VpAr. Using the interior L°° bounds on p^r, Ujv (Propo- 
sition Hj3 points 2 and 5), J \ Vpn\ 2 is uniformly Lipschitz in compact sets of ]0, T[. Hence we can 
conclude that 

pT p rp N-k N . 

(45) lim / / \Vp N (t,x)\ 2 dtdx - — V / | Vp N (t i7 x)\ 2 dx = 0. 

v • x Jo JTrf M rr* Jr d 

(Remember that since r n \ 0, we also have k^/N ~ j 0.) We can conclude that 
limF A r(pAr, Ujv) = \im F N (p, v) = F(p, v ) as N — ► oo. 

We show now that this implies that (pn, v n) and (p, v) are close to each other. For this we use 
the dual formulation of the problem (Proposition 18 .5|) . This formulation needs an admissible solution 
(p, J = pv,p) satisfying (fTTH ITT ] IT2" | fi"3 |) . We take here (p,v,p) the optimal solution of Problem 11.31 
which is admissible for Problem 18.11 Using the fact that Fpt(p,v) is close to Fn(pn,vn) for iV large, 
for any e, 5 > there exists N, i/) N , q e N (with N = N(S)) such that 

N-k 

Y\ / \V PN (U)\ 2 dx + 5 

< e + 5+ dpd t tl) e N + dpv -Vif>jf 
Jd 

rp N-k 

+ mJ2 p(U)qN{U)+Vp(U)-W N {U)-\WN(U)\ 2 /2dx. 



1 /" T 

< -J Pn\vn\ 2 dtdx + — 



i=k 

This eventually yields 

N-k 



f 1 T" p ^ 

j D 2^ ~ V ^' 2 dtdX + 2N^ dX J T d 2 |Vp(tl) ~ Vg ^ )|2 rfx 



< e + 5. 



For fixed iV , (Vip e N , q f N ) being a maximizing sequence for the dual problem, it will converge to 
(vn,Pn) as e — > (see (jSJ)), therefore we obtain 



\ f p\v- v N \ 2 dtdx + — Y] [ \Vp(ti) - Vpn(U)\ 2 dx < 5, 

1 J D ~T k JT d 

and therefore this goes to as iV — > oo. Using the same procedure we can also get that 



Pn\v — vn\ 2 + |VpN — Vp| 2 dtdx — > as N — > oo. 



Now we show that the product Pat^at converges to pt>: using the equicontinuity property of the 
sequence ||pat(£, OHz^ m [ r > ^ — r ] fo r any r < T/2 and any 1 < k < oo, (see Proposition 18. 6|) . the 
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sequence pn converges strongly in L k ([r,T — t] x T d ) for any 1 < k < oo. Moreover remember that 
from Theorem EH v G L 2 ([r, T - r] x T d ). Then 

/ \pv - p N v N \ 

J[T,T-r]xT d 

< / Pn\v n - v\ + \v\ \p N - p\ -> 0. 

•/ [t,T— t] X T d 

Hence, pjvfAr converges strongly to pv in L} oc (D). Using the uniform L^ C (D) bound on p^, vn (points 
2 and 5 of Proposition 18. 6j) . we get that PnVn converges strongly to pv in Lf oc (]0, T[xT d ) for any 
1 < p < oo. 

Proof of Theorem 12.31 The theorem is obtained passing to the limit in the Proposition 18.61 The 
point 2, 3, 4, 5 remain true when we let N go to oo. The other points will be shown at the end of 
the paper. 

□ 

Then from Theorem l2~3l if p ,Pr & Ld+?(Y d ) we have p G L°°([0, T]; (T d )) and doing as in 
Lemma 13.0.11 Vp G L°°([0, T]; L 2 (T d )). This bound shows that the Riemann sum 
If ^1^=1 fjd IVp(^> a; )| 2 ^ a; converges to J D \ Vp\ 2 dtdx and this allows us to take r = in the previous 
proof and to conclude the proof of Proposition 18.71 

□ 



8.3 Regularity properties of the time discretized solution: Proof of 
Proposition 18.61 

In this part N is fixed and for sake of simplicity we drop the subscript N. We consider 
(p = Ap + 1, v — V0) solution of the Problem 1)8.1)1 . 

8.3.1 Preliminary: Construction of a special solution 

First we begin to show the consistency with the optimal transport Problem 18.31 let U < s,t < t i+1 , 
we denote 

(46) * S)t (x) = (t-s)<f)(s,x) + \x\ 2 /2. 

The function $ Sit goes from M. d to M if we extend to a periodic function on all of M. d . Note also 
that for any p G Z d , V<fr Sj t(. +p) — V +p. If s = U and t = t.- l+ i we denote (resp. 
instead of $ Sit (resp. $t, s ). Note that is discontinuous at times U and thus 

$ M+1 (x) = \x\ 2 /2 + ^(tt,x), 
T 

&i,i-i(x) = \x\ 2 /2 - —<j>(t i} x). 

We also have 

(47) v(s, x) = V0(s, x) = (V* a , t (x) - x), 

t — s 

which is well defined on R d /Z d . 
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In the first lemma, we will see that V«fr Sj t pushes forward p(s) on p(t) minimizing the cost 
L d p(s, x)\m(x) — x\ 2 dx among all m pushing forward p(s) on p{t) (that we denote hereafter 
m # p(s) = p{t)) and that <fr s t coincides with its convex hull dp(s) a.e.. 

Then in the second lemma we will show that we can consider a solution for which every $ s t is 
convex. This point that may seem to be a direct consequence of optimal transport (the fact that 
the optimal transport is given by the gradient of a convex function) needs from our point of view 
a careful proof: indeed we only know that the gradient of <& st will coincide dp(s) a.e. with the 
gradient of a convex function, but since the optimality equation links 3>i_i,i, and p(ti) it must 

be checked that can consistently be taken convex. The convexity will then allow us to consider 

the second derivative of since a convex function is almost everywhere twice differentiable, and 

then to make rigorous the inequality (J4*0|) and its consequences. 

First for any / : R d — > R we denote /* its Legendre transform defined by 

(48) f*(y) = sup{yx-f(x)}. 



The convex hull of / is (/*)* (or /** in short). It is equivalently defined as the supremum of all 
convex functions smaller than /. We will show the following lemma: 

Lemma 8.8 Let U < s,t < U+i and <& t>s , $> st be defined as above. Then 

<& s t > <&l s with equality dp(s) a.e., 
$ fl , t = $** dp(s) a.e., 
V^ st = V$:* dp(s) a.e., 
V*3 # p(a)=p(f). 

Proof. We know that cf) is the limit of a smooth sequence <ft e satisfying the the constraint 

N 



T rc| N 

i=l 



Between ti and t i+ i we have 



(49) <9 t e + ^|V0 e | 2 <O. 

y — ^ 

Thus if t > s considering 7(0") = x + (a — s) and using (}4T?|) we find 

t — s 

A[0 e((7)7((7 ))] = ^ (a;7((T)) + p^.V0 e (a,7(a)) 
da t — s 

< d,W,lW)) + \ + |V^(ff,7(ff))| 2 

< lj» 



2 V t- s 



2 t 



Integrating from s to t we find 



\y — *-H 

<Pe{t,y) < inf{0 e (s,x) + — -}, 



12 

2(t-s) 



12 



\y — ^ 

(s,x) > sup{(p £ {t,y) - — -}. 

y s ) 
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Letting e go to 0, it follows from (jUaj) that 

(so) & t ,s(y) > • y - ®sA x )} = KM, 

X 

(51) > 

This is the first point of the lemma. The crucial point is the following: for dp(t) a.e. y we have 

I 1 2 

y) = inf {(j)(s, x) + 2 _ }, 

or equivalently 

$t, s (y) = sup{x • y - $ s ,t{x)} dp(t) a.e.. 



Indeed take a smooth sequence (4> e ,p e ) such that 

N 

X 



1 T N 



i=l 

that maximizes the dual problem, i.e. such that 

T f 

-/T d 
A/"—] 

+ ^7^/ PPe(U,x) + Vj> ■ Vp e (ti,x) - -\\7p e \ 2 (U,x) dx 

> K N -e 2 . 

In view of (jUj) , being a maximizing sequence of the dual problem implies the following: 

r i i - 

/ P\dt<t>e + ^<Pe\ 2 + ■jr^2S t=t .p e (t i )\ dtdx -f as e -> 0, 
^ D ~ »=l 

which in turn implies that 

/ pd t 4>e + pV0 ■ V0 e dt'dx 

pit, x)(f> t (t, x) — p(s, x)<f> e (s, x) dx 

1 ft r i 

p|V0| 2 dt'dx = W2(p(s),p(t)) as e -> 0. 

where the first line comes from the mass conservation equation (j2J) satisfied by the pair(p, v = V0) 
and the last identity comes from the fact that between t and s the problem coincides with the optimal 
transport Problem 18.31 

If m is a mapping realizing the optimal transport of p(s) onto p(t) then m#p(s) = p(t) implies 

p(t,x)(j) e (t,x) — p(s,x)(p e (s,x) dx — / p(s, x)((f) e (t,m(x)) — (f) e (s, x)) dx. 
jd Jjd 
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Using that 

\x ?/| 2 

<f> e (t, x) - <f) e (s, y) < — - _ ^ 
and that from the optimality of m we have 



2{t- 

we obtain by taking the limit e — > that 



If 1 

r / p(s,x)\x - m(x)| 2 dx = W£(p(s),p(t)), 

- s) J T d t-s 



\x — m(a;V 2 



which is equivalent to 



(f)(t,m(x)) = <p(s,x) H _ ^ dp(s) a.e., 



<f>{t,y) = <fr(s,x) + - 2 (f-s)^ dp ^ a ' e " 



We remind the reader that m is invertible dp(t) a.e. and m -1 can be defined as the (dp(t) a.e. 
unique) mapping realizing the optimal transport of pit) onto p(s). This implies that if t > s, we 
have: 



\y — x 

(pit, x) = inf{— r + 0(s, y)} dp(t) a.e., 

y 2{t — s) 

®t,s = (*,,*)* dp{t) a.e., 

the two lines being equivalent through equation (}46|) . As a supremum of affine functions, any Legendre 
transform is convex. Here, <& t ,s coincides with a convex function dp(t) a.e. and is above this function 
dx a.e. from ([50)1 . Since $** is the convex hull of s it follows that 

$** = $ M dp{t) a.e., 

from which it can be deduced that 

/ (|x| 2 /2 x) dx+ [ i\y\ 2 /2~^liy))pit,y) dy 



= (t — s) p(s,x)(j)(s,x) — p(t,x)(j)(t,x) dx 
= WHpis),pit)). 
This implies that V<fr** t# p(s) = p(t) (see [S]). Note also that if we set 

x) = -L [(| ■ | 2 /2 + (t - s)<f>(s, •))** (*) - |x| 2 /2] 

we obtain that 

/ p(t,s)|V0- V(p\ 2 it,x)dx = 
for a.e. t. The proof of Lemma f8. 81 is complete. 



□ 
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We are going to use the previous lemma to construct a new sequence of solutions for which the 
potentials 3?t,i+i are convex. This will allow us to define dp(ti) a.e. the second derivative of $^+1. 



rp A-! 



This special solution will turn out to be the viscosity solution of d t ip + ^IV^I 2 + -j-j $t=uP = 



Remember that from Proposition 18 A\ <f)(tf) satisfies 

T 

(52) (p(tj, x) - 4>(ti , x) < -jrP{U, x) dx a.e., 

T 

(53) <P{tt,x) - </>(ti ,x) = -—p{U,x) dp(U) a.e.. 



i=l 



Consider the new solution ip defined by 



(54) ^ t = 0,x) = ^ 



l - + ^(t = 0,-)) (x)-\x\ 2 /2 



(55) on }ti, t i+ i [, tp{t, x) = inf <^ — — — + tp{tf, y) 

y I2{t-ti) 



N 

(56) ip( t +,x) = - 



Exchanging px and po, (p(T — t),p(T — t), —<p(T — t)) will be the corresponding optimal solution, 
and we introduce also ip constructed by the same procedure from (p(T — t), —<p(T — t)). 

Lemma 8.9 

1. For almost every t e [0, T], ip(t, .) coincides with <ft(t, .) dp(t) a.e. and (p,v = Vif),p) is solution 
of Problem \8.1\ 

2. Mi G [0..N — 1], (t,x) — > ip(t,x) and (t,x) —ip(t i+ i + U — t,x) are both viscosity solutions 
(and subsolutions) of d t ip + llVV'l 2 = on [ti,t i+ i], 

rp n-i 

3. ip is the viscosity solution of d t ip + llV^I 2 + -j-j ^t=uP = on [0, T] in the sense of / T^9j) . 

i=i 

4- finally —ip(T — t) = ip{t) holds dp(t) a.e.. 
Proof. We denote 

(57) ¥.,t(x) = \ x \ 2 / 2 + (*- S M S ' x ) for s ' * e fc-i, *<+i]> 

and ^1,1+1, as well. Let us now prove by induction the following: 

For all 1 < i < N — 1, ip(t~,x) > <p(t~,x) with equality dp(ti) a.e.. 

Equation fl^U) implies that * ,i = $0*1 • Tlien (J23) implies that *i i0 = *o,i = *o*T = $ o,i - *i,o 
with equality dp(ti) a.e. from Lemma IHtHl The equality $qT = $01 comes from the fact that $J 2 
is convex as a Legendre transform and that the the Legendre transform is an involution on convex 
functions. Thus from (J57j) we get that ip(t±) > (f)(t^) with equality dp(ti) a.e.. This proves the 
property for i — 1. 
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Suppose now that ip{t i ) > ) with equality dp(ti) a.e.. Then we have, using (jozjl 

|x| 2 /2 + -V(t-,x)-— p(ti,x) > |x| 2 /2 + x)-— p(t h x) 

> |x| 2 /2 + ^0(t+ x), 

where both inequalities are equalities dp(U) a.e., hence ^i+i > <£**+i from (}56|) . We now show that 

= <&** +1 dp(ti) a.e.. 

Indeed from Lemma ()8.8)1 . we know that = dp(U) a.e., therefore the convex function 

is below |x| 2 /2 + ,x) — j^p(ti,x), with equality dp(tj) a.e.. We conclude since ^i+i is 

pinched between |x| 2 /2 + j^ip(t^,x) — ^pit^x) and that coincide dp(tj) a.e.. 



Then 



\J> • , , • = \J>* , < <b**.* , = „ < * • , , ■ 



where we have used that the Legendre transform is decreasing with respect to functions, is an 
involution on convex functions and from Lemma 18.81 for the last inequality. This implies ip(t~ +1 ) > 

We have 

* M+1 = dp(ti)a.e., 
V¥ M+ i = V$** +1 dp{ti) a.e., 
V* M+ i#p(^) = V$** +1# p(ti) = p(t i+l ). 

This implies that ^j+i = (^*i+i)* dp{ti+i) a - e - : f° r this use the identity 

c(x) + c*(Vc(z)) = x ■ Vc(x) 

that holds for any Lipschitz convex function c. This will also imply ( = (<fr£* t )* dp(t) a.e. for all 
t E [ti,ti+i]. Hence we have proved that 4>(t^ +1 ) = ip{t~ +l ) dp(ti + i) a.e.. In view of Lemma 18781 we 
also have for all t G [0,T], ip(t) = <p(t) dp(t) a.e.. This proves the first point of the lemma. 

On ]tj,tj+i[ we have d t ip + ^\Vip\ 2 = satisfied in the viscosity sense from the Hopf-Lax formula 
(j33J). To see that —ip{t i+ i + U — t,.) is also a viscosity solution to this equation, first note that 
= |£| 2 /2 + jfip(ti,x) is convex. Then we have 

\x — y\ 

tp{t,x) = m{{tfj(U,y) + — — — } 
v 2{t — ti) 

■ r r i / \ \ x ~~ y\ 2 \ x ~~ y\ 2 ^ 

= mt{ip{ti, y) H — + 7— -} 

v 2{t i+1 -ti) 2(t i+1 -t) 

for some 7 > 0. Now observe that a(y) = i/j(ti,y) + )^~ y [, ; is convex as well as (3(y) = 7 Jf 
and thus we can apply the Fenchel-Rockafellar duality Theorem that says that if a, (3 are convex 
continuous we have 

mi{a(y) + f3(y)} = sup{— a*(x) — (3*(—x)}. 
y x 
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Computing a* gives -. j and we can check that 

W+l I ' 



3/ — y\ 

if;(t,x) = sup{ip{ti +1 ,y) - — -}, 

which says exactly that t —* — ip{t i+ i + i, — t) is the viscosity solution of d t if) + || Vip\ 2 = on 
Then we check that 



if)(tf +l ,x) 



y [(i • r/2 + •) - 7 ^f 1 r(x) - ixi 2 /2] 



< ^(* i+ i,a;) 



Tp(U,x) 



with equality pjv(^+i) a e., which shows that ^ is a subsolution of 



(58) 



1 T JV_1 

^ + -|V^| 2 + -^5^ = 



To see that if) is the viscosity solution of (f58]). just notice that our definition of t/> is the following: 

n 2 /2 - = (i-r/2+^(*+,or 

= d-i a /2+^(*r,o-Sp(^or, 



thus 



and we obtain 



T 

TO+1> X ) = illf W< > - xrW. ?/) + 
2/ iV 



F - 1/1 



2— 



mf ^(*r.7(*r))-+ 



/ -^^<W(^7(^)) + ^ItIV) 
i=i 



da 



where T is the set of all continuous paths with j(ti +p ) = x, and more generally that 
(59) ^(t,x) = 

■ r j, N-l 1 



inf <^(s,7(s)) + 

j(t)=x 



i=l 



da > . 



This defines the viscosity solution of (}58|) . and proves the point 3 of the lemma. 

The proof of the point 4 follows simply from the point 1: ip(t) and <p(t) coincide dp(t) a.e., and if> 
is built from —<p(T — t) in the same way as if) is built from 0, therefore if>(t) and —<f){T — £) coincide 
dp(T — t) a.e., and we conclude. This achieves the proof of Lemma 1531 

□ 

Remark. The time reversibility property is not valid for any viscosity solution. Actually this is 
true before occurrence of shocks. This is what one says when we decompose 0(0, y) + ^t-t- ) as ^ e 
sum of two convex functions: this means that one can continue the rays further without developing 
shocks. Thus we see that our variational solution does not develop shocks in the interior of the time 
interval. 
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8.3.2 Proof of the bound on A<fi 
The function ip satisfies 

^ + 1|V^| 2 + ^W(^)<0, 
d tP + V • (pVVO = 0. 

From now we consider that := ip an d thus is convex for any t{ < s,t < t i+ \. We are going to 
prove the following lemma: 

Lemma 8.10 The functions are C 1 ' 1 at every density point of p(ti) for 1 < i < N — 1, 

with C 1,1 norm bounded by C(N, d). Moreover there exists a set Si of full measure for p(ti) such that 
everywhere in £i, (j)(tf) and 4>(t~) are twice differentiable and the following holds 



(60) 



A0(t+ x) - A(f>{tr, x ) < ^(1 -pfax)). 



Proof. Using (jSSl E2J) we get that 



(61) 
(62) 



T 2 



&i,i+i(x) + & iti _i(x) = \x\ - —p(ti,x) dpiti) a.e. 

T 2 

&i,i+i(x) + &i ti -.i(x) < \x\ 2 - j^p(ti, x) dx a.e.. 



Thus, dp(ti) almost everywhere, the convex function 3>i,j+i + is tangent from below to 

| ■ | 2 — j^p(U, ■). The Poisson equation satisfied by p{U) being only true in the distribution sense we 

need to introduce a finite difference version of the Laplacian: 

Lemma 8.11 Let Cd be the volume of the unit ball ofM. d and bd be the d — 1 dimensional Hausdorff 
measure of the unit sphere. Let p be continuous and 



AM = % 



1 



1 



c d h d 



p{y)dy - p(x) 

8B(x,h) 

p(y)dy-p(x) 



B(x,h) 



with the constants kd, Id chosen so that both operators converge to the Laplacian for smooth functions 
as h — > 0. Then 

1. if Ap > C in T d in the distribution sense then A^p > C and A* h p > C everywhere in T d , 

2. if Ap G L 2 (T d ) then, up to extraction of a subsequence in h, A* h p — > Ap dx almost everywhere, 
finally, point 2 still holds even if p is not continuous. 

Proof of claim 1. We take / solution of 

/ _ P\dB(x,h) = 0, 

Af = C < Ap. 

Observe that the boundary condition has a meaning since p is continuous. Then from the max- 
imum principle f(x) > p(x) in B(x,h), and thus Ahf < AhP since / = p on dB(x,h). But 
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f(x) = C\x\ 2 /2d + g with Ag = 0; using the fact that the average on a sphere of an harmonic 
function equals its value at the center of the sphere, we get that AyJ = C. To obtain the inequality 
for A* h just integrate over h. 

Proof of claim 2. If suffices to show that A* h p converges strongly in L 1 to Ap if Ap G L 2 (T d ). 
The Taylor formula gives for dx almost every x: 



A* h p(x) 



1 l d 

c d h d h 2 



[ [ (1 - a)y* ■ [D 2 p{x + ay) - D 2 p{x)} -ydody 

J\y\<h Jct=0 



+ [ \v l ■ D 2 p(x) ■ ydy 

J\y\<h 1 



= E 1 (x) + E 2 (x). 

Then £2(2?) is equal dx almost everywhere to Ap(x) and Si converges strongly to in L 2 since 
p G W 2 ' 2 from elliptic regularity. 

□ 

Hence, since Ap(tj) = p(ti) — 1 G L 2 (Y d ) (from Proposition 18. 4|) . up to extraction of a subsequence 
in h, for almost every x G T d , A* h p(ti, x) converges to p(ti, x) — 1. Applying A£ to (|UT1 at a point 
where equality (jfiTjl holds we get that 

(63) A*$ i>i+1 + A* h ^ x <2d- ^A* hP (U). 

Therefore we have dp(U) a.e. 

limsupA*$ M+1 + A*^_x < 2d+ — (1 -p(U)). 

A convex function, say /, such that A* h f(x) is bounded as h goes to is C 1,1 at x, since the trace 
controls the norm of a positive matrix. Moreover the left hand side is nonnegative, implying that 
dp{ti) a.e. 

2dN 2 

p{U,x) < 1 + -j^-. 

This yields a L°° bound on p (that depends on N), and implies from elliptic regularity that p(ti) is 
(ji,a f or ever y < a < 1. Since convex functions are continuous and p(U) is continuous, equality 
(|HT|l holds now on a closed set of full measure for p{ti) and thus at every density point of p(U). Now 
we can apply point 1 of Lemma IS. Ill to obtain that A* h $> iji+ i + A^j^x is bounded at every density 
point of p(ti), implying the expected C 1 ' 1 bound. Convex functions are twice differentiable almost 
everywhere, and thus almost everywhere A* h $> iji+ i — ► A$j ii+1 and A^ i)4 _j — > A$i fi _i. We define Si 
to be the set of all points x such that 

1. x is a point of Lebesgue differentiability for p(tj) where p(ti,x) > 0, 

2. A* h p(ti, x) — > p(ti, x) — 1 as h —> 0, 

3. and are twice differentiable, 

4. equality (joTj) holds. 
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The set Si is of full measure for and by removing sets of measure to Si one can impose the 
condition V^j+i^) = S i+ i. Then at every point of Si we have 



□ 



(64) 0<A^ i>i+1 (x)+A^ 1 (x)<2d+—(l-p(t i ,x)), 

2dN T 

(65) -— < A0(t+ x) - Acf>(tr, x) < -(1 - p(t,, x)). 

This implies the following bound for p, <p: 

\\p{U)\\L~<C{d)N\ 

for every x G S u ||0(ti, x)|| c i,i < C{d)N/T. 
The proof of Lemma 18.101 is complete. 



Construction of the characteristics Remember that is given by: 

& i>i+ i{x) = |x| 2 /2 + ^0(t+,x), 

and that for s G [ti,t; + i] (resp. s G , & ti ,s is given by 

*uA x ) = \x\ 2 /2 + (s-U)<J)(tt,x), 
( resp. & ti ,s( x ) = \ x \ 2 / 2 + ( s -ti)4>(ti,x)). 

We know that in Si, $^+1, are convex and twice differentiable. Thus for any s G]£j_i, U+i[, 3?t i)S 

is twice differentiable and D 2 <& ti ,s is invertible. This implies (see the Appendix on convex functions 
of [21]) that D 2 <& st . exists at point V$ tiiS (x) for x G Thus we can define £ s = V$ iijS (^), 
this definition makes sense pointwise in £ i7 and pointwise in S s , <p s is twice differentiable. Note 
also that from Lemma 18.91 the definition S s = V<&t ijS (£j) and S s = V<&t i+liS (£j+i) are consistent 
since V<&i,i+i(£i) = Then we can define a trajectory (x s ,s G [0,T]) as follows: starting from 

X 8Q G £ So for any s g]0,T[, we define x s = V<f> So ,s(xo) for any s in the same interval [tj,tj + i] as 
So, and we proceed similarly in other intervals. Thus we define a flow S(s,t, x) that gives at time 
s the position of the particle located in x at time t. This flow H(s, t, 2) is defined everywhere on 
£t, and S(s,t, £ t ) = S s . We may denote (x s , s G [0,T]) a trajectory and it will be understood that 
Xt G St, Vt G]0, T[. 



Conclusion of the proof Now we bound A(f> along a trajectory: for (t, s) G [tj,tj+i], we have 

(t-s) (A<f>(t)(x t ) - A(f>(s)(x s )) 
= (t-s) (A0(t, V* S)t (z.)) - A0(s, x,)) 
= 2rf-A$ Sji (x)-A$ M (V$ M (x s )), 

but this is negative since we have the relation 

D 2 $> s ,{x) = [D 2 ^ s }-\V^ t (x)), 

and therefore 

2d - A$ s , t {x) - A$ M (V$ M (x s )) = 2d- ^(Aj + 1/AO < 0, 
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where the Aj are the eigenvalues of D 2 $> Stt (x s ) well defined on £ s . Thus we conclude first that for 
every x s e £ s and x t = V<& Sjt (a: s ) 

(66) (A0(t, x t ) - A0(s, ac a )) • (t - s) < for t 4 < s, t < t i+1 . 

Then we obtain a quantitative estimate of the decay of A0, between ti and t i+ i\ we take s = t^t = t i+ i 
in the previous inequality, from the convexity of x — > x + l/x we have 

d d 

y2(Xj + 1/Xj) > d(A/d + d/A) where A = ^ Xj = A$ M+1 

3=1 3=1 

l -{A-df + 2d 



A 

;gA0(t+ Xi )y +2d 



since j;A(j)(tf) + d = A$ iii+ i. Hence 



A" 

N jfA<f>(tr,Xi) + d 

Using ((HQ)) we obtain 



iV ^ d + ±A(/)(tt,x. 

We know from fH that A$ iii+1 < 2d + T 2 /iV 2 thus A0(t+,x;) < f (d + T 2 /iV 2 ). It follows that 
(2 + jfA(f)(tf, Xi) < 2d + ^ < 3d for TV large enough. We finally obtain the following bound for N 
large enough: 

A0(t+ +1) x m )<A0(t+ x^ + Kl lAo «t-'-i»- 



N V 3d 

This is a discrete version of the differential inequality < 1 — g^@ 2 and we will conclude that 

A(p(U) < C(d)(l + -) in £ for any 1 < i < N - 1. 

ti 

This will be a consequence of the following lemma: 

Lemma 8.12 Let (-X" n )ne[o.jv] be a sequence defined by X 6 R and sttc/i t/iat 

^n + l<X n + ^(l-X 2 /C 2 ). 

T/ien for n > 1, 

x <c+° 2 t - Tn 
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Proof. The proof is by induction. For n — 1, we look for 



rnax{X + -(l-X 2 /C 2 )}. 

This yields 

Xi < T/N + NC 2 /{AT) <C + C 2 /t n = C + C 2 N/{nT) 

for N > T/C, 1 < n < 4. Then suppose that n > 4, X n < C + C 2 /t n . The function 
x -> x + ^(1 - x 2 /C 2 ) is increasing for a; < NC 2 /(2T). For iV > AT/C and n > 4, we have 

NC 2 /{2T) >C + C 2 /t n = C + NC 2 /{nT). 

Therefore we have (for N > AT/C, n > 4) 

X n+l < C + C 2 /t n + ^(l-(C + C 2 /t n ) 2 /C 2 ) 
= C + C 2 /t n (l-2T/(CN)-T/(Nt n )). 

In order to conclude we need 1 — 2T/(CN) — T/(Nt n ) < n/(n + 1). But we have 
1 — T/(Nt n ) = 1 — 1/n < n/(n + 1), therefore we conclude the proof of the lemma. 

□ 

Using then the transformation <f)(t,x) — > —<p(—t,x), that transforms the solution of Problem 18.11 
in another solution of 18. II exchanging p and px, we get that Acj) N (t i ) > —C(d)(l + tJ-t) in Si for 
any 1 < i < N — 1. Since we know from ()66)1 that t — > Acf)(t,Xt) is decreasing between and ti + i 
we can conclude that there exists for each t a set of full measure for dp(t) on which <p(t, .) is twice 
differentiate, and where the following equality holds: 

(67) - C(d)(l + ^— f ) < A<j> N (t,x) < C(d)(l + - t ). 

The first part of Proposition 18.61 is proved. 



8.3.3 Proof of the Lf^QO, T[xT d ) bound on p 

We begin by writing the Monge- Ampere equation that links p(tj) to p(t i+ i) 

T 

p(t i+2 , x i+2 ) det(J + —D 2 (p(tf +1 , x i+1 )) = p(t i+1 , x i+1 ), 
T 

p(U, x^ det(J - —D 2 <p(tr +11 x i+1 )) = p(t i+1 ,x i+1 ). 

This equation makes sense since, from Lemma 18.101 at x s , <ft is twice differentiable, and we use 
[2U Proposition A2]. 

Now using the domination of the geometric mean by the arithmetic mean we have 

det(J + ^£>W + i)) < (1 + ^A0(t+ 



hence 

(68) < (i + ^mt^ 



p(t i+1 ,X i+ i) < n T 
p(t i+2 ,Xi +2 ) ~ dN' 



(69) ^!<^i± 1 )< (1 __L A , (im , I , +l)) , 
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We deduce first the following: 



1 p(t i+ i,x i+ i) T ss d 

< 72 r~ - U ~ T77 A 0(^+i ; ^+i)) • 



dN 

Note that we also have 



(1 + -§A0(t+ Xi )) d ~ p(ti, m) ~ dN 



<™) 3Tg^ ^^ (1 -^^ 



for ti < s,t < t i+ i, x s G £ s , and this implies using (|57|) that along a trajectory s — > x s , log(p(s, x s )) 
is Lipschitz: for all t 1? t 2 £ [ r ? ^ — r], 

(71) | Mp^J - log(p(t 2 , s t2 )| < C(l + — - - — -)\t 2 - t x \. 

r(T - t) 

Taking the logarithm of (J68{ l69|) we obtain 

log(p(t i+2 , X i+2 )) + l0g(p(tj, Xi)) - 2 log(p(t i+ i, Xj+i)) 

> -dlog(l + J^A0(t+ 1; x m )) - dlog(l - ^A<f>{tr +1 ,x i+1 )) 



T 

N' 
T 2 



(72) > j^(p(t i+l ,x i+1 )-l), 

where at the third line we have used the concavity of the log and at the last line we have used (|6U|): 

T 

A(j)(tf +1 , x) - A0(t i+1 , ar) < — (1 - dp(t i+1 ) a.e.. 

We fix r g]0, T/2[. For any trajectory x s ,s G]0,T[ with x s G £ s for all s, log(p(s, x s )) is uniformly 
Lipschitz with respect to s in [r, T — t] from (jTTj) and log(p(s, x s )) remains finite in [r, T — r]. 
Moreover (J?2"J) holds at every time £j. Using this we claim an unconditional bound for p(s,x s ) for 
t < s < T — T. 

Proof of claim: The sequence (©(ti) = log p(ij, £i))i<i<jv-i satisfies a discretization of the differential 
inequality 

(73) 0>exp0-l. 
From (|71|l . we have the a priori bound 

|e(t i+ i)-e(tj)| <c(t)t/n 

for in [r, T — r]. We argue by contradiction: take < r < T/4 and suppose Q(U ) > M, with 
t io G [r, T/2]. Then choose M so large that M - C(r)(T - It) > M/2. Thus, on [r, T — r] we have 
© > M/2 from the a priori bound above. This implies that 

[©(* i+ i) - ©(ti)] - [Q(ti) ~ e(Vi)] > T 2 /iV 2 (exp(M/2) - 1) 

on [t, T — r]. We then choose M large enough so that Q[(t-t)n]/t ~ ®[(t-t)n]/t-i > C(r)T/N in 
contradiction with the a priori bound above. 

Applying the same argument after having changed t in T — t gives the bound on [T/2, T — r]. 
Note that this proof does not depend on the initial and final values of 0. 
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8.3.4 Time continuity of p 

Remember that in Theorem 12.21 we have proved that p G C(]0, T[; L p (T d )) for any p G [1, §[• We 
have now an unconditional bound on p in L^ c (]0, T[xT d ). Thus the strong time continuity in every 
L p (T d ), 1 < p < oo follows and the point 2 of Proposition 18. 61 is proved. 

8.3.5 Lipschitz bound for log -)Ili fc (Trf)) 

Since is twice differentiate at x s G £ s , we can use the identity 

pit, x t ) det(7 + (t- s)D 2 4>(s, x s )) = p(s, x s ), 



and it implies 



We have thus 



— f [p(t,x)] k dx = -(k-l) [ [p(t,x)] k A(j)(t,x) dx. 
at Jjd Jjd 



j t \W,-)\y < ^\\A(f)(t r )\\ L ^ dp(t)) \\p(t,-)\\ Lk , 



hence using point 1 of Proposition 18.61 we get that 



-CW 1 + t) ^ I lo s Ollzw) < c(d)(i + 1 



t' ~ dt V 3 „l^x-„ - ■ 

This proves the point 3 of Proposition 18.61 

8.3.6 Displacement convexity of functionals of p 

Here we show that given pw(t, x) solution of Problem 18. 1[ the functions f Jd pjv(£, x) \og(piy(t, x)) dx, 
Jjd[PN(t, x)] k dx, k G [l,+oo[ converge to convex functions of t G [0, T]. We drop the subscript N. 
Let 4 denote the convective derivative d t . + V0 ■ V. and -C = (^) 2 - The density p satisfies (|72*jl 
which is the finite difference version of 

d 2 

(74) — log p(t,x) >p(t,x)-l. 



Note that we do not include all terms, since one could show that our solutions satisfies 

+ p(t, x) - 1. 



d 2 . . , 1 ' 



^ 2 logp(t,x)> d 



^p(M) 



The term we have omitted here allows to extend the displacement convexity of (k — 1) /[p(£, x)] k dx 
down to A; > 1 — l/d in the non-interacting case (cf [21]). However in the gravitating case, the 
additional term will limit this property to positive values of k — 1. What could probably be proved in 
this case is that between k = 1 and k — I — 1/d the functional (k — 1) J[p(t,x)] k dx is semi-convex. 
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Formal proof We first give a formal proof, eluding the fact that we only have a discrete version 
of the differential inequalities. The lemma given below will justify the calculations. Using (J71j) with 
the identities 

^ X) = ^t~x) ' ^ X) ' 2 + ^ X) ^ ^ ^ ' X) ' 
^[p(t,x)] fc = fc(* - l)[p(t,x)] fc - 2 ||p(t,x)| 2 + fclXt.x)]*" 1 ^^), 
we obtain for k > 

^[p(t,x)] fc > fc 2 [p(t,x)] fc - 2 |^p(t,x)| 2 + k[p(t,x)]\p(t,x) - 1). 



Noticing that for all smooth F we have 
dt 2 



d 2 

d p(t,x) — (F(p(t,x))) dx, 



p(t, x)F(p(t, x)) dx 

and applying this to F(p(t,x)) = [p(t,x)] k , k > 0, we get 

/ [p(t,x)] k+1 dx> [ k([p(t,x)) k+2 -[p(t,x)] k+1 )dx>0 
at jfd j^d 

Indeed, using Jensen's inequality twice we have 



[p{t,x)] k+1 dx > / [p{t,x)] k dx 



k+1 
k 



k 



[p(t,x)} k dx) / [p{t,x)} k dx 
Jd J Jjd 



> 1 / \p(t,x)]* dx 

for k > 1 to conclude. 

This convexity property combined with the unconditional bound for p in L°°([r, T — r] x T d ) yields 
a uniform Lipschitz bound for \\p(t, -)||L fc i n [ T iT — r] for any 1 < k < oo. (For the case k = +oo 
this is because of (|7T|) ). 

Rigorous proof 

Lemma 8.13 Let fc 6e greater than 1. Then p, the limit of the sequence p n solutions of Problem \8.1l 

satisfies 

d 2 r r 



, [p(t,x)] fe dx>(k-l) I [p(t,x)] k+1 - [p(t,x)] k dx 

in the sense ofV'(0,T). 

Proof. We consider p solution of Problem IS. II We find after some elementary calculations that 

Pi+i(x i+ i) + Pti(^-i) - 2p k (xi) 
> kp k i ~ l (x i )(p i+ i(x i+1 ) + pi_i(xi_i) - 2pi(xi)) 

+ -k(k - l)p k ~ 2 (xi) ((p i+1 (x i+1 ) - Pi(xi)) 2 + (pi_i(xi_i) - Pi(xi)) 2 ) 
+ 0{N- 3 ). 



51 



The term 0(N~ 3 ) depends on || A0 (t) 1 1 x,<» (d/3(t)) » ||p(*)IU°° an d is therefore uniformly bounded in com- 
pact sets of ]0,T[. Then we use (|72j). to get also by simple calculations (using the concavity of the 
logarithm) 

p i+ i{x i+1 ) + Pi-i(xj_i) - 2pi(xi) > -j-ppi{xi){pi{xi) - 1) 

1 / (p i+ i(x i+1 ) - pi(xj)) 2 | (pi-i(^-i) - pi{xj)f 



2 V Pi( X i) Pii X i) 

Combining these inequalities, we get 

pIi(x w ) + ftiM - 2 P k M) > ^-M^))ip^) ~ !) + k2 Q + o(n- 3 ), 

with Q a positive quantity The term 0(N~ 3 ) is uniform on every interval [r, T — r], for r > 0. We 
integrate over T d , this yields 

/ (pffi + pK-apT 1 ) rfx>^ jk(pt 2 - P n dx + o(N- 3 ). 

We now choose a test function <p e C£°(]0,T[). Denoting Z}^ the second incremental quotient 
(Dlf(x) = h~ 2 (f(x + h) + f(x -h)- 2f(x)), with h = T/N), we have 

rp N-l / r \ rr N ~ 1 r 

nY,[ d I \p&> x )f = jiT, bte> x t D lv&) 
i=i ^ ^ ' i=i ^ 

-> / [p(t,x)]V(i) dtcfe 
in 

which proves the lemma. 

□ 

The convexity of J p log p follows the same lines, and we skip the proof. 
8.3.7 Proof of the W 1 >°°{[t,T - r] x T d ) bound for <p 

Hereafter we use again the subscript N for the solution of Problem 18. II while V0,p is the solution of 
Problem 11.31 

From Lemma HUE <Pn can be given by the Hopf-Lax formula (J59|) that we recall here: 



r ,t[ T N-i 1 

(75) <f> N (t,x)= inf <<Ms,7(s))+/ -— V 5 a = t .p N (a, 7(a)) + -|7lV 

7(*)=a: I Js iV ~~[ 1 



do 



where we still use U = jj- This formula is valid for any < s < t < T. We are going to prove the 
following lemma that will yield the result when letting N go to +00. 

Lemma 8.14 Let (p^ be defined by \7dj) . If\VpN(Ui ")l — ^( r ) f or an V^i £ [ r >^~ r ] with I (t) < +00 
/or r > 0, and p^(t^-) has mean value 0, then 

|V<Mt,x)| <C(r), 

|^(t, *) - 0at(s, s) I < C(r) (^(1 + E( iV(t ~ s) )) + |t - g| 



for any < r < t, s < T — r, dx a.e. x G T d , where E(-) denotes the integer part. 
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Proof. Let 7 be a minimizer in the infimum (|75)1 . Let 7(0") = 7(0") + i~{z ~ x )- Then 



<f>N(t,z) < (f> N (s,i(s)) + 



2V-1 



N 



r ^5 CT=ti ^(a, 7 (a)) + ^| 2 (a) 
i=i 

+ / -^^^=^(^7M) + d7lV) 

" /s L i=l 



da 



da 



+ / W")ljf^ + ^$*'+(i + ^(^))'WI*-»l- 

We choose here s, t such that t — s > r/2, s > r/2 and t <T — t. We will have 7 uniformly bounded 
in L 2 ([s,t]) by a constant C(r). The second line is equal to <j)^(t,x), the third line is bounded by 
C(r)^r an d we g e ^ that 



||V0Ar|| L o O ([ T)T _ T ] xT d) < C(r). 
Now we choose s,t G [r, T — r]. By taking 7(0") 



x 



<j>N(t,x) < 0at(s,x) + C(t)jj(1 + E( ^l ^ )). Using now that V0tv is bounded by C(r) we get 



in (|75|) . we obtain that 



4>N(t,x) - (f) N (s,x) 

r t N(t — s) \y - x\ 2 

> inf <^ (j) N (s,y) - (j) N (s,x) - — (1 + E( ))|bjv|Uo 0([s t]xTd) + — 

> inf (-^(1 + EC^^))!^^!!^^^^,) - CCr)^ - a;| + 1 ?/ " ' 



yei d 



N 



T 



2(t - s) 



> -C'(r)[^(l + E(^L^l)) + \t-s\ 



This proves the lemma. □ 
We can use this lemma to conclude the W 1,00 ([t,T — t] x T d ) bound for <fi. Indeed, we know 
already that \p N (t,x)\ < C(r) for t G [r,T-r]. Since Ap N = p N - 1, we thus have p N (t) G W 2 > k (T d ) 
for 1 < k < 00, uniformly on [r, T — r]. The Sobolev embedding Theorem then yields 

lbiv|U^([T,T--r];C' 1 (T d )) — C( r )- 



We let then N go to +00 so that jjE( N ^ s ' ) — > t — s. We already know that pn converges to p in 
C(]0, T[; L p (T d )) for p < 00. From the bounds obtained we will also have p^ converging uniformly 
to p on compact sets of ]0,T[xT d . Finally, extracting a subsequence if necessary, <p N will converge 
to uniformly on compact sets of ]0, T[xT d . 



8.3.8 Convergence to viscosity solutions 

For a given smooth 7 : [s, t] — > T d compute 



<f N {t,x) = (f) N (s,j(s)) 



^(t,x)=^s, 1 (s)) + 



J s i=\ 



da, 



-p( ( T,7(a)) + -| 7 | 2 (a) 



da. 
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Combining the mass conservation equation (J2J) and the Poisson equation, we get 

d t p N = A-\-V- ( PN V<j) N )). 
From the Lf£ c (]Q, T[xT d ) bounds on p N , V(pN and elliptic regularity, we get that for any d < k < oo 

d tPN (t, .) e W^ k (T d ) c C a (T d ) 



a 



uniformly for t G [r, T — r]. We also already know that Pw(t) is uniformly Lipschitz in space for 
t G [t,T — t]. Hence we have <j>Jf(t, x) — > 7 (t, x) since (4>n,Pn) converge uniformly to (0,p) in every 
compact set of ]0,T[xT d and from the bound just obtained on p. Since we can choose 7 to remain 
bounded in L 2 ([s,t]) and thus 7 bounded in Ca([s, i]) we conclude that 



inf {^(f,*)}- inf {^(t,*)} 

ll7lli,2 [Sit] <C,7(t)=a! ll7lli2 [Sit] <C,7W=a! 



as — > 00. Hence we get that 
(76) <j>(t, x) = inf 

j(t)=x 



(s,7(s)) + 



-p(<T,7(«r)) + -|7| 2 (a) 







(icr 



and is the viscosity solution of d t <t> + \ \ V0| 2 + p = on every [s, t] c]0, T[. 

Remark. In particular, <p is a subsolution, therefore admissible for the time continuous dual 
problem. Hence, if 4>n maximizes the time discretized dual problem, it converges to the maximizer 
of the time continuous dual problem, uniformly on compact sets of ]0,T[xT a! . 

8.3.9 Reversibility and C 1,LlogL regularity for cf) 

We can define the backward solution ip starting from 4>(T) by ip{T) = (f>(T) and 



(77) 



ip{s,x) 



sup 

7(s)=x 



-p(o-,7(o-)) + -|7| 2 (^) 



da 



for any < s < t < T . This amounts to take for ip the limit of the sequence — iPn{T — t) in 
Lemma 18.91 Therefore t — ► —ip(T — t, •) is the viscosity solution of d t (p + i|V<^| 2 + q = with 



g(t, ■) = p(T — t, •). Then from Lemma f8. 91 we know that ip = <p dp a.e.. Since t/> and are Lipschitz 
continuous, the set <fi = ip is a closed set of full measure for p(t) for each t G]0,T[. This means that 
on the set that contains the dynamics (i.e. the support of p), the solution of our Hamilton- Jacobi 
equation is reversible, a property which is generically false for viscosity solutions. This property is 
close to several regularity properties obtained in [Tlj on the Aubry set in weak KAM theory. This 
implies regularity on the support of p: 

Using that in the time interval [s,t], p — 1 + Ap is uniformly bounded, we have from elliptic 
regularity (see for example) Vp uniformly log-Lipschitz, (that we will denote by p G C 1,LXogL ) 
i.e. V(x,y) G T d with \x - y\ < ~,V£ G [r,T-r], 



\Vp(t,x)-Vp(t,y)\<C(r)\x-y\log( 



\x 



Then consider an optimal path 7 in (|76|). Perturb it in 7(0") 



l{°) + l^t(z 



X) 



We have 



(j>(t,z)<Tl x (z) = <t>(t,x) + / 7(a) 



.r 



p(a,7 



a — s 
t- s 



\z 


— x 2 


\t 


-S| 2 



c/cr 



-p(o-,7) 



do - . 
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The second integral is C 1,LlogL with respect to z, while the first is a quadratic polynomial in z. 
Moreover we have Tf x (z) > <ft(t, z) and T£ x (x) = 0(£, x), which means that T^ x is tangent from above 
to 4>{t, •) at x. 

In a similar way, we use (|77j) to obtain at every point (x,t) a C 1,LlogL function T xt tangent from 
below to if). The log-Lipschitz constants are uniformly bounded in [r, T — r]. 

If follows that at every point of the set {if> = 0}, (or ip) is pinched between two functions whose 
gradients are log-Lipschitz. This implies the following, which is the point 7 of Theorem 12.31 

Proposition 8.15 For every t G [r, T — r], r > 0, t/iere exists a closed set St of full measure for 
pit) such that is differentiate at every point of St, and 

V{x,y)eS t ,\x-y\ < 1/2, | V0(t, x) - V<f>{t,y)\ < C{r) \x - y\ log{-^—). 

\x - y\ 

Proof. We define the set S t = {x G T d , (f>(t, x) = ip(t,x)}. We choose two points (x, y) G S t . 
Using that <fi is pinched at x between two C 1,Ll ° sL functions, we have at x, for all z G T d , 

\<f>(t, z) - <f>(t, x) - V(f)(t, x)-iz-x)\< C\x - z\ 2 log( , 1 , ), 

\x — z\ 

and we have the same changing x in y. Combining the two inequalities, we get that 
\(V<j)(t,y)-V<f)(t,x))-(z-x)\ 
< C (\x - z\ 2 \ogi-. — - — r ) + \y - 3| 2 log(-| — - — r ) + \y — x\ 2 \ogi-. — - — r ) ) . 

V f — z \ \v ~ A \v ~ x \ J 

Taking the maximum value of the first line among all z such that {\z — x\ = \y — x\} (note that this 
implies \y — z\ < 2\y — x\), we get that 

| V0(t, y) - V0(t, x)\ < C'\x- y\ ]og(—L-). 

\x-y\ 

This proves the proposition. 

□ 

8.3.10 Estimates up to the boundary 

If p T G L k (T d ) with k > d, then p remains Lipschitz in space up to t = T, hence one can choose 
(f) G W 1 ' 00 ^, T] x T d ). Instead of (JZHJ) one can have defined by (JUJ), hence if p G L k (T d ) } one can 
choose G W 1 ' oo ([0,T - r] x T d ). If p and p T are both in L k (T d ), with k > d then, using (JZZj) one 
can choose first such that 0(0, •) G W 1 ' 00 ^). Then we obtain that G W liOO ([0,T] X T d ) using 
the following result that can be found in ^2] : 

Proposition 8.16 Let be the viscosity solution on [0, T] x T d of 

^0+^|V0| 2 + p = O, 
0(t = O,-) = 00, 

withpe L oo {[0,T};W 1 ' oo iT d )) and O G lU 1 ' 00 ^). T/ien G ^'^([O,^ x T d ). 
This last result ends the proof of Theorem 12.31 

□ 
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